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Abstract 


The  Pawnee  Grassland  Eddy  Correlation  Dry  Deposition  Project  is 
described.  Instrumentation,  methods  of  analysis,  and  initial  data  and 
research  findings  are  presented.  Data  from  this  eddy  correlation 
system  show  agreement  with  (a)  previous  observations  of  deposition 
velocities  for  atmospheric  ozone,  NOz,  and  NOx;  (b)  micrometeoro- 
logical  theory;  and  (c)  micrometeorological  site  characteristics. 
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MANAGEMENT  IMPLICATIONS 

The  deposition  of  sulphur  (S)  and  nitrogen  (N)  com- 
pounds is  the  focus  of  a  large  research  program  being 
carried  on  currently  throughout  the  United  States.  The 
objective  of  this  research  program  is  to  determine  what 
loadings  of  S  and  N  (in  kilograms/hectare/year]  will  cause 
adverse  ecological  effects  and  to  quantify  these  effects. 
Furthermore,  research  on  forested  ecosystems  suggest 
that  ozone,  a  pollutant  photochemically  generated  from 
emissions  of  nitrogen  oxides  and  hydrocarbons,  must 
also  be  included  with  S  and  N  in  determining  forest 
health. 

Depending  upon  climate  and  the  ambient  atmospheric 
environment,  S,  N,  and  ozone  may  be  deposited  by  either 
dry  "processes"  or  by  precipitation.  This  paper  discusses 
only  the  dry  deposition  process. 

At  first  glance  measuring  dry  deposition  might  seem 
to  be  an  easy  task;  for  example,  simply  put  out  a  bucket 
and  measure  what  falls  in  when  it  is  not  raining.  How- 
ever, dry  deposition  measurements  are  not  simple 
because  the  surfaces  on  which  dry  deposition  occurs 
may  themselves  actively  be  removing  the  pollutant. 
Therefore,  simple  passive  monitoring  is  not  adequate  for 
a  full  description  of  the  dry  deposition  process.  Some 
understanding  of  the  nature  of  the  underlying  surface 
is  also  required. 

The  most  accurate  techniques  for  measuring  dry 
deposition  employ  aerodynamic  methods.  Essentially 
these  involve  measuring  how  much  pollution  passes  by 
a  fixed  level  in  the  atmosphere.  This  procedure  works 
best  when  the  vegetation  is  uniform  and  the  terrain  is 
flat.  The  research  reported  on  in  this  paper  represents 
a  test  of  the  aerodynamic  method  known  as  eddy  cor- 
relation under  circumstances  where  it  is  expected  to 
work  well. 

Eventually,  we  plan  to  use  eddy  correlation  to  make 
dry  deposition  measurements  in  typical  alpine  and  sub- 
alpine  ecosystems.  These  eddy  correlation  measure- 
ments will  then  be  used  to  confirm  simpler  technologies 
for  routine  measurements  of  the  levels  of  S,  N,  and  ozone 
present  and  depositing  to  forested  lands. 

INTRODUCTION 

Awareness  that  acid  precipitation  and  dry  deposition 
of  acid  gases  and  aerosols  can  have  potentially  signifi- 
cant impacts  on  the  biota  at  the  earth's  surface  has  grown 
steadily  in  recent  years  (NAPAP  1987,  Wetstone  1980). 
However,  the  processes  by  which  pollutants  are  trans- 


ported to  and  subsequently  interact  with  the  biota  and 
the  soils,  and  the  influence  that  these  pollutants  can  have 
upon  the  biota,  remain  largely  unknown.  In  an  effort  to 
better  understand  these  processes  the  USDA  Forest  Serv- 
ice, Rocky  Mountain  Forest  and  Range  Experiment  Sta- 
tion (USFS/RMFRES),  Atmospheric  Deposition  Research 
Project,  has  undertaken  a  dry  deposition  monitoring  pro- 
gram. This  paper  describes  that  program,  outlines  the 
data  analysis  methods  used,  and  presents  some  initial 
results  from  our  measurements. 

This  report  specifically  details  efforts  to  implement  the 
eddy  correlation  method  of  estimating  the  dry  deposi- 
tion fluxes  of  ozone,  N02,  and  NOx.  This  micrometeor- 
ologically  based  method  uses  a  direct  measurement  of 
the  covariance  between  the  fluctuating  vertical  velocity 
and  the  fluctuating  concentration  of  the  acid  gas  to  deter- 
mine the  vertical  flux  (or  rate  of  deposition  per  unit  of 
surface  area)  of  the  pollutant  to  the  underlying  surface. 
The  eddy  correlation  technique  is  not  the  only  method 
available  for  estimating  fluxes  [Wesely  and  Hicks  (1977) 
and  Businger  (1986)  give  reviews  of  this  and  other 
meteorological  methods];  but,  it  does  have  several  advan- 
tages over  the  other  approaches.  It  is  the  most  fundamen- 
tal of  the  different  meteorological  methods  because  it 
requires  the  fewest  assumptions  and,  at  present,  is  also 
the  only  method  capable  of  providing  a  direct  measure 
of  the  vertical  flux  of  a  quantity  to  or  from  the  underly- 
ing surface.  It  is  important  to  note  that  the  eddy  correla- 
tion flux  data,  although  taken  at  a  fixed  point,  can  be 
considered  to  be  the  integration  of  the  surface  fluxes 
from  some  distance  upstream.  A  disadvantage  of  the 
eddy  correlation  method  is  that  it  is  the  most  demanding 
technically  because  dry  deposition  rates  are  particularly 
difficult  to  measure.  For  example,  Sehmel  (1984)  lists  51 
separate  variables  that  can  affect  dry  deposition  rates. 

The  paper  first  gives  the  physical  and  biological 
description  of  our  research  site  and  an  overview  of  the 
instrumentation  used  before  addressing  the  errors  asso- 
ciated with  measurement  systems  in  general  and  this 
system  in  particular.  Then,  initial  results  of  the  data 
analysis  are  presented.  Also  included  are  several  appen- 
dixes that  detail  important  aspects  of  our  overall  effort. 

PAWNEE  RESEARCH  SITE  DESCRIPTION 

The  Pawnee  Project  eddy  correlation  research  site  is 
50  km  northeast  of  Fort  Collins,  CO,  at  the  Central  Plains 
Experimental  Range  (CPER),  a  research  area  adjacent  to 
and  west  of  the  Pawnee  National  Grasslands.  The  CPER 
is  a  6,500-ha  area  of  shortgrass  prairie  that  is  ad- 
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ministered  by  the  Agricultural  Research  Service  of  the 
USDA.  Also  located  at  the  CPER  is  a  long-term  ecological 
research  (LTER)  shortgrass  steppe  research  station 
(Halfpenny  and  Ingraham  1984).  The  LTER  station  is 
managed  by  Colorado  State  University  (CSU)  and  is 
funded  by  the  National  Science  Foundation.  The  exact 
CPER  location  of  the  eddy  correlation  research  site  is 
within  a  6-ha  grazing  exclosure  in  Weld  County  (section 
27,  1,641  m  elevation,  latitude  40°28'23"N,  and  longitude 
104°45'15"W.)  (Fig.  1).  This  exclosure  was  established 
in  1969  as  an  appropriate  and  "representative"  site  for 
meteorological  instrumentation.  The  area  is  level, 
bordered  on  all  sides  by  gently  undulating  plains,  and 
is  not  near  major  ridges  or  bluffs  that  would  seriously 
influence  measurements. 

Vegetation  on  the  CPER  is  dominated  by  shortgrasses 
(64%),  succulents  (20%),  and  half-shrubs  (8%).  The 
immediate  area  surrounding  the  eddy  correlation  instru- 
mentation is  dominated  by  blue  grama  grass  {Bouteloua 
gracilis),  the  major  plant  community  type  of  the  short- 
grass  steppe  area  in  eastern  Colorado  (Lauenroth  and 
Milchunas  1985).  Blue  grama,  a  warm  season  grass  that 
exhibits  the  dicarboxylic  acid  biochemical  pathway  of 
C02  fixation,  has  an  optimum  photosynthetic  temper- 
ature (when  soil  water  is  not  limiting)  of  between  26°C 
and  33 °C  (Brown  and  Trlica  1974).  Also  present,  but  less 
common,  are  Bouteloua  gracilis-Agropyron  smithii  (blue 
grama-western  wheatgrass)  plant  communities.  Western 
wheatgrass  is  a  "cool  season  grass"  (greatest  growth  dur- 
ing the  spring  and  early  summer)  that  is  infrequent  in 
the  research  vicinity.  However,  Stipa  comata  (needleand- 
thread),  also  a  cool  season  grass,  is  common  and  nearly 
a  codominant  in  the  immediate  vicinity  of  the  Pawnee 
Project  instrumentation  tower.  The  most  frequent  forbs 
in  the  area  are  Opuntia  polyacantha  (pricklypear  cactus), 
Gutierrezia  sarothrae  (snakeweed),  and  Artemisia  frigida 
(fringed  sage).  Table  1  provides  further  details  on  species 
composition  within  the  exclosure  site. 

The  maximum  total  green  leaf  area  index  (LAI)  at- 
tained from  studies  of  various  vegetation  types  at  the 


Table  1.— Plant  species  in  the  immediate  area  of  the  NOAA/eddy 
correlation  experimentation.3 
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Figure  1.— Site  location  and  vertical  cross  section  of  the  instrument 
shelter,  chemical  instruments,  and  meteorological  tower. 


Cactacae  Opuntia  polyacantha  Haw. 

Chenopodiaceae  Chenopodium  album  L. 

Chenopodiaceae  Kochia  scoparia  (L.)  Schrad. 

Commeliniaceae  Tradescantia  occidentalis  (Britt.)  Smyth 

Compositae  Artemisia  frigida  Wi lid.* 

Compositae  Chrysothamnus  nauseosus  (Pall.)  Britt. 

Compositae  Cirsium  undulatum  (Nutt.)  Spreng. 

Compositae  Grindelia  squarrosa  (Pursh)  Dunal 

Compositae  Gutierrezia  sarotharae  (Pursh)*  Britt.  +  Rusby 

Compositae  Senecio  tridenticulatus  Rydb. 

Compositae  Thelesperma  filifolium  (Hook.)  A.  Gray 

Compositae  Tragopogon  dubius  Scop. 

Convolvulaceae  Evolvulus  nuttallianus  R.  +  S. 

Cruciferae  Lepidium  densiflorum  Schrad. 

Cyperaceae  Carex  eleocharis  Bailey 

Cyperaceae  Carex  filifolia  Nutt. 

Graminae  Aristida  longiseta  Steud.* 

Graminae  Bouteloua  gracilis  (H.B.K.)  Steud.** 

Graminae  Buchloe  dactyloides  (Nutt.)  Engelm.^ 

Graminae  Sporobolus  cryptandrus  (Torr.)  Gray* 

Graminae  Stipa  comata  Trin.  +  Rupr.* 

Malvaveae  Sphaeralcea  coccinea  (Pursh)  Rydb.* 

Polygonaceae  Eriogonum  effusum  Nutt. 

^Nomenclature  taken  from  Flora  of  the  Great  Plains,  1986,  T.M. 
Barkley  etal.  (editors),  Great  Plains  Flora  Association,  University 
Press  of  Kansas,  Lawrence,  Kansas,  1392  p. 

bFrequency  of  occurrence  indicated  as  follows: 

'Common  in  tower  vicinity. 

"Very  abundant  in  tower  vicinity. 

All  others  occur  infrequently. 

CPER  was  0.55  (Knight  1972).  Peak  LAI  was  reached 
about  mid-June  and  from  60%  to  80%  of  total  LAI  was 
blue  grama.  In  1971  the  greatest  LAI  increment  occurred 
between  April  20  and  June  15,  with  a  steady  decline  after 
June  15.  This  seasonal  pattern  of  total  LAI  production 
for  steppe  vegetation  was  similar  to  the  pattern  of  total 
amounts  of  chlorophyll  (Rauzi  and  Dobrenz  1970).  Pat- 
terns similar  to  these  can  be  expected  for  any  year,  but 
the  rate,  time,  and  amount  of  LAI  or  chlorophyll  produc- 
tion for  a  given  year  will  reflect  different  precipitation 
and  vegetation  patterns.  For  example,  the  LAI  produc- 
tion in  the  vicinity  of  the  Pawnee  Project  meteorological 
tower  may  have  a  bimodal  peak  because  of  the  frequent 
needleandthread  in  addition  to  the  dominant  blue  grama. 

Standard  meteorological  data  has  been  gathered  by 
CSU  within  the  6-ha  grazing  exclosure  for  several  years. 
In  addition,  the  exclosure  includes  a  8-m2  weighing  soil 
lysimeter.  Seventy  percent  of  the  311  mm  of  mean  an- 
nual precipitation  occurs  between  April  and  September. 
This  low  rainfall,  generally  flat  to  gently  rolling  terrain, 
the  absence  of  local  air  pollution  sources,  along  with  the 
long-term  meteorological  data  records,  make  it  a  good 
site  for  attempting  to  measure  dry  deposition. 

The  exclosure  site  described  above  is  also  a  pollutant 
deposition  research  facility  hosting  three  separate  acid 
deposition  monitoring  programs: 

1.  National  Atmospheric  Deposition  Program/National 
Trends  Network  (NADP/NTN).  Wet  and  dry  deposition 
chemistry  measurements  are  collected  and  analyzed  at 
a  central  NADP  laboratory  in  Champaign,  IL.  Data  from 


2 


this  site  (NADP  site  no.  06222)  are  available  from  March 
1978  to  present. 

2.  NOAA/ATDD  Dry  Deposition  Inferential  Monitor- 
ing (DDIM).  This  site  is  no.  11  of  the  NOAA/ATDD  core 
dry  deposition  sites.  Filter-pac  chemical  analyses  of  site 
samples  collected  during  the  same  sampling  period  and 
frequency  as  the  NADP/NTN  wet  chemistry  samples  are 
conducted  by  the  USGS  laboratory  in  Atlanta,  GA.  Dry 
deposition  chemistry  is  coupled  with  micrometeor- 
ological  data  to  infer  dry  acid  species  fluxes  (Hicks  et 
al.  1985].  Dry  deposition  chemistry  and  associated 
meteorological  data  are  available  from  July  1986  to 
present. 

3.  USFS/RMFRES,  Atmospheric  Deposition  Research 
Project,  Pawnee  Eddy  Correlation  Acid  Deposition 
Research  Site.  Ambient  air  concentrations  and  eddy  flux 
measurements  of  N02,  NOx>  ozone,  heat,  and  momen- 
tum have  been  collected  during  two  short-term  data  col- 
lection periods  from  mid-December  1986  through 
January  1987  and  from  mid-May  1987  through  July  1987. 
This  last  program  is  the  subject  of  this  paper. 

Our  eddy  correlation  system  requires  simultaneous 
measurements  of  the  fluctuations  in  the  vertical  wind 
speed,  heat,  and  chemical  concentrations.  These  meas- 
urements are  made  at  a  rate  between  10  and  20  Hz. 
Hence,  both  the  wind  sensor  and  chemical  sensors  must 
have  relatively  short  response  times  (see  appendixes  A 
and  B  for  a  discussion  of  response  times).  Meeting  these 
fast  response  time  requirements  requires  very  sensitive 
meteorological  sensors  and  chemical  analyzers.  The 
chemical  analyzers  for  this  measurement  program  were 
provided  by  the  University  of  Denver  (appendix  A). 
These  instruments  are  housed  in  an  instrument  shelter, 
and  ambient  air  was  brought  to  them  using  a  high- 


volume  sample  pump  and  a  large-bore  Teflon  tube.  A 
more  detailed  description  of  these  sensors  is  given  in  ap- 
pendix A  and  by  Ray  et  al.  (1986),  Wendel  et  al.  (1983), 
and  Wesely  et  al.  (1982).  Table  2  lists  the  sensing  instru- 
ments and  analyzers  integrated  into  the  overall  eddy 
correlation  program. 

During  the  two  field  periods,  all  eddy  correlation 
measurements  were  made  using  a  3-axis  Gill  propeller 
anemometer  for  wind  velocity  and  velocity  fluctuations, 
and  a  fast  response  platinum  resistance  thermometer  for 
temperature  fluctuations.  The  Gill  anemometer  was  cor- 
rected in  real-time  for  non-cosine  response  (Massman 
and  Zeller  1988).  During  one  short  test  period  in  Decem- 
ber 1986,  a  one-dimensional  sonic  heat  flux  sensor  was 
used  to  compare  heat  fluxes  acquired  with  the  Gill  ane- 
mometer and  the  fast  response  platinum  resistance  ther- 
mometer. We  also  plan  to  replace  the  Gill  anemometer 
in  future  experiments  with  a  3-axis  sonic  anemometer. 

Data  collection  and  processing  were  accomplished 
using  a  personal  computer  interfaced  with  an  analog-to- 
digital  conversion  board.  Eight  channels  are  available 
on  the  board  for  data  acquisition  and  all  eight  channels 
(table  3)  were  sampled  at  ca.  14  Hz.  Output  measures 
include  means,  variances,  covariances,  and  other  micro- 
meteorological  parameters  as  listed  in  table  3.  The  soft- 
ware used  for  real-time  data  processing  was  adapted 
from  McMillen's  (1986)  basic  flux  program  for  eddy 
correlation  in  nonsimple  terrain.  The  compiled  basic  pro- 
gram is  capable  of  running  the  entire  eddy  correlation 
measurement  system  unattended  for  several  days.  Tasks 
accomplished  in  real-time  by  this  program  include  real- 
time data  acquisition  and  processing,  raw  data  storage 
if  required,  statistical  summaries  at  the  end  of  the  as- 
signed averaging  period  (usually  1/2  hour),  instrumen- 


Table  2. 


-Instrumentation  used  at  Pawnee  Acid  Deposition  Research  Site  for  eddy  correlation 

measurements.3 


Instrument  or 
model 

Typo 

Parameter 

Response 
time 

Range 

NOx  analyzer 

chemiluminescence 

NOx 

0.9  s 

0.07  ppb  -  1000  ppm 

CAAM  2 

chemi  luminescence 

N02 

0.7  s 

0.1-1000  ppb 

CAAM  1 

chemiluminescence 

o3 

0.7  s 

0.2-400  ppb 

TECO  49 

UV  adsorption 

o3 

2  min 

1-500  ppb 

R.M.  Young  UVW 

Gill  propeller 
anemometer 

wind  velocities 

(b) 

0.3-30  m/s 

AIR,  Inc. 

platinum 
resistance 

temperature 

(b) 

-50-  +  50°C 

Atmos.  Tech.  Inc. 

3-axis  sonic 

wind  velocities 

(b) 

0-20  m/s 

Campbell 
Scientific 

vertical  sonic 
(CA27T) 

heat  flux 

±16°C/±4  m/s 

AT&T  6300  + 

PC  computer 
(80286) 

data  collection 
and  processing 

Data  Translations,  Inc. 

DT2805  board 

A/D  conversion 

aThe  use  of  trade  and  company  names  is  for  the  benefit  of  the  reader;  such  use  does  not 
constitute  an  official  endorsement  or  approval  of  any  service  or  product  by  the  U.S.  Depart- 
ment of  Agriculture  to  the  exclusion  of  others  that  may  be  suitable. 

bSee  appendix  B. 
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Table  3.— Data  acquisition  and  measured  parameters  at  Pawnee  Acid  Deposition  Research  Site. 


Channel 

Parameter 

Outnut  fpvprv  1/2  hr\ 

vuipui  fcvci  y   lit  illy 

1 

w:  vertical  wind  speed 

w~,  W2 

2 

u:  horizontal  wind  speed  component 

IT,  u72 

3 

v:  horizontal  wind  speed  component 

v~,  V2 

4 

T:  temperature 

T,  T'2 

5 

03  (Teco  49):  ozone 

03  ,  03'2 

6 

03  (CAAM1):  ozone 

03  ,  03'2 

7 

N02  (CAAM2):  nitrogen  oxide 

N02  ,  N02'2 

8 

NOx  (nitrogen  oxide) 

NOx  ,  NOx'2 

Momentum  flux 

w'u',  w'v' 

Heat  flux 

w'T' 

Ozone  flux 

w'03' 

N02  flux 

w'N02' 

NOx  flux 

w'NOx' 

Horizontal  wind  speed 

(u2  +  V2)1'2 

Wind  direction 

arctan(v/u) 

* 

Friction  velocity 

(-w'u')1'2 

Drag  coefficient 

(-w'u'  )lu2 

Horizontal  wind  rotation  angle 

V 

Vertical  wind  rotation  angle 

e 

Sample  frequency 

#  of  samples/sample  period 

N 

'Derived  from  data  recorded  on  first  8  channels. 


all  others,  quite  unlike  what  was  expected  from  the 
1972-75  data  set. 

EDDY  CORRELATION  SYSTEM  ERRORS 

Our  dry  deposition  eddy  correlation  system  employs 
a  variety  of  instruments  and  support  equipment;  still,  no 
data  system  is  without  sources  of  error.  For  example, 
Businger  (1986)  reviews  and  comments  on  several  recent 
attempts  by  various  researchers  to  measure  dry  deposi- 
tion of  gases  and  particles.  He  notes  that  up  to  the  pres- 
ent there  has  been  no  discipline  among  the  various 
researchers  to  ensure  that  their  data  have  all  the  ingre- 
dients for  an  accurate  determination  of  dry  deposition, 
and  he  presents  10  specific  sources  of  error  that  he 
claims  are  found  in  most  field  programs: 

1.  Failure  to  account  for  the  effect  of  heat  and  water 
vapor  fluxes; 

2.  Sampling  errors  due  to  insufficient  averaging  time; 

3.  Inadequate  sensor  response  characteristics; 

4.  Sensors  separated  by  too  great  a  distance; 

5.  Random  noise  not  accounted  for; 

6.  Errors  inherent  in  the  horizontal  uniformity  and 
steady-state  assumptions  due  to  entrainment, 
advection,  and  nonstationarity; 

7.  Errors  due  to  wrong  sampling  height; 

8.  Inadequate  fetch; 

9.  Errors  related  to  deliquescences  of  particles; 
10.  Errors  due  to  flow  distortion  around  sampling 

probes  and  mounting  structures. 
In  order  to  make  routine  measurements  of  dry  deposi- 
tion fluxes,  it  is  necessary  to  examine  these  and  other 
errors  and,  where  possible,  account  for  them  in  the 
resulting  data.  We  begin  by  examining  errors  that  can 


tation  zero  each  half  hour  to  ensure  accurate  chemical 
measurements,  and  instrument  span  calibration  each 
half  hour  (not  yet  physically  implemented). 

The  physical  layout  of  the  eddy  correlation  site  is 
shown  in  figure  1.  The  relative  directional  location  of 
the  instrument  shelter  site  to  instrument  tower  was  based 
on  meteorological  data  collected  during  1972-75,  which 
showed  the  210°  azumuth  as  one  of  the  least  likely  wind 
directions  (Peterson  and  Parton,  undated).  However,  the 
analysis  of  one  year  of  10-m-height  wind  data  at  the 
Pawnee  Project  tower  (fig.  2)  clearly  shows  that  during 
daylight  hours  southerly  wind  directions  dominate  over 
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Figure  2.— Pawnee  site:  a,  annual  wind  rose;  b,  daytime,  growing 
season  wind  rose. 
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arise  from  sensor  limitations  in  size  and  response,  and 
instrument  separation  and  mismatch. 


Flux  Loss  Corrections 

Methods  similar  to  those  of  Moore  (1986),  as  detailed 
in  appendix  B,  were  used  to  estimate  the  flux  losses  with 
the  Gill  propellers  in  the  sensor  array  and  then  the  data 
were  corrected  for  these  errors.  The  heat,  momentum, 
and  chemical  flux  losses  for  stable  and  unstable  at- 
mospheric conditions  as  a  function  of  wind  speed  are 
given  in  figures  3,  4,  and  5,  respectively.  The  solid  lines 
assume  matched  sensors  and  the  dashed  lines  on  these 
figures  assume  mismatched  sensors  (if  no  dashed  line 
appears  then  sensor  mismatch  was  not  significant).  By 
matched  sensors  we  mean  that  the  two  separate  types 
of  instruments  used  to  make  the  flux  measurements  have 
identical  response  characteristics  and  mismatched 
sensors  mean  they  do  not.  The  separation  distances  used 
for  these  computations  were  36.0  cm  between  the  ver- 
tical Gill  propeller  and  the  temperature  sensor,  63.0  cm 
between  the  vertical  Gill  propeller  and  the  horizontal  Gill 
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Figure  3.— Estimates  of  heat  flux  loss  as  a  function  of  wind  speed 
and  atmospheric  stability  associated  with  using  the  Gill  propeller 
anemometer  for  measurements  of  the  vertical  velocity  fluctua- 
tions. Solid  lines  assume  matched  sensors,  dashed  lines  assume 
mismatched  sensors. 
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Figure  4.— Estimates  of  momentum  flux  loss  as  a  function  of  wind 
speed  and  atmospheric  stability  associated  with  using  3-axis  Gill 
propeller  anemometer  system  for  measuring  the  fluctuating 
velocity  components. 
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Figure  5.— Estimates  of  chemical  flux  loss  as  a  function  of  wind 
speed  and  atmospheric  stability  associated  with  using  the  Gill 
propeller  anemometer  for  measuring  the  vertical  velocity  fluctua- 
tions. Solid  lines  assume  matched  sensors;  dashed  lines  assume 
mismatched  sensors. 

propellers,  and  37.8  cm  between  the  vertical  Gill  pro- 
peller and  the  mouth  of  the  intake  tube  leading  to  the 
chemical  sensors.  For  matched  sensors,  this  system  will 
underestimate  the  heat  flux  by  at  least  20%  to  30%  for 
unstable  conditions  and  about  50%  for  stable  conditions. 
Likewise,  the  momentum  flux  is  underestimated  between 
about  10%  and  60%,  depending  upon  wind  speed  and 
stability.  These  latter  results  are  quite  consistent  with 
Horst's  (1973b)  estimates  of  the  momentum  flux  loss  that 
can  be  expected  with  Gill  propellers.  Finally,  the  ozone, 
N02,  and  NOx  flux  losses  are  between  30%  and  90%, 
again  depending  upon  stability  and  wind  speed.  The 
NOx  flux  loss  is  slightly  higher  than  that  of  the  other  two 
sensors  because  of  its  slightly  slower  response  time. 
These  large  flux  losses  result  from  the  slow  response 
times  for  both  the  Gill  propeller  and  the  chemical 
sensors. 

The  dotted  lines  on  figures  3  and  5  give  the  estimates 
of  the  heat  flux  losses  and  the  chemical  flux  losses  caused 
by  mismatching.  For  the  heat  flux,  mismatching  can  in- 
crease the  flux  loss  by  another  8%  for  an  unstable  at- 
mosphere and  another  15%  for  a  stable  atmosphere.  For 
the  chemical  fluxes,  at  low  wind  speeds  the  estimates 
of  flux  losses  are  also  increased.  For  a  stable  atmosphere 
the  flux  loss  is  an  additional  8%  for  the  chemical  fluxes, 
and  for  an  unstable  atmosphere  it  is  about  another  5%. 
As  figure  5  also  shows,  mismatching  has  little  effect  on 
the  chemical  fluxes  at  moderate  wind  speeds,  because 
the  vertical  Gill  propeller  has  a  response  time  com- 
parable to  that  of  the  chemical  sensors.  Furthermore,  this 
figure  also  shows  that  as  the  wind  speed  increases, 
mismatching  again  begins  to  contribute  additional  errors 
because  of  the  reduction  of  the  response  time  of  the 
vertical  Gill  propeller  at  high  wind  speeds  relative  to  the 
response  time  of  the  chemical  sensors.  Similar  compari- 
sons were  made  for  the  momentum  flux  losses.  However, 
the  response  times  of  the  horizontal  and  vertical  Gill 
propellers  were  not  different  enough  to  introduce  any 
significant  mismatching  errors. 

In  addition  to  the  above  errors  associated  with  the  Gill 
UVW,  two  other  matters  relating  to  its  use  for  estimating 
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turbulent  fluxes  and  variances  are  relevant.  The  first  in- 
volves the  need  to  correct  the  turbulence  data  gathered 
with  these  instruments  for  their  lack  of  a  perfect  cosine 
response.  Massman  and  Zeller  (1988)  discuss  this  at 
length  and  provide  an  algorithm  that  can  be  used  to  make 
this  correction  in  real-time.  They  show  that  without  this 
correction,  turbulent  heat  and  momentum  fluxes  are 
underestimated  by  about  30%. 

Unfortunately,  the  second  problem  cannot  be  fully  cor- 
rected. Because  these  propellers  tend  to  stall  whenever 
the  flow  is  greater  than  about  88°  off  axis,  they  will 
always  underestimate  turbulent  fluxes.  This  is  most 
serious  with  the  vertical  propeller,  which  senses  a  field 
that  is  always  fluctuating  through  the  propeller's  stall 
region.  This  problem  will  be  most  severe  whenever  the 
horizontal  velocity  is  significantly  greater  than  the  ver- 
tical velocity.  Finkelstein  et  al.  (1986)  discuss  in  greater 
detail  the  limitations  of  these  sensors  for  turbulence 
measurements. 

Again,  the  methods  outlined  in  appendix  B  were  used 
to  estimate  the  flux  loss  errors  associated  with  using  the 
3-axis  sonic  anemometer  in  place  of  the  Gill  anemometer. 
For  this  case,  the  separation  distances  are  6.35  cm, 
0.0  cm,  and  30  cm  between  the  center  of  the  vertical 
sonic's  acoustic  path  and  the  temperature  sensor,  the 
center  of  the  acoustic  path  of  the  horizontal  sonic,  and 
the  mouth  of  the  intake  tube  to  the  chemical  sensors, 
respectively.  The  solid  curves  and  dashed  curves  of 
figures  6,  7,  and  8  give  the  results  of  the  computations 
for  the  flux  losses  for  the  sonic  anemometer  correspond- 
ing to  figures  3,  4,  and  5,  respectively.  In  general,  the 
3-axis  sonic  anemometer  provides  much  better  estimates 
of  the  fluxes  of  heat  and  momentum  at  all  wind  speeds 
and  the  fluxes  of  the  three  chemical  species  at  low  wind 
speeds  than  does  the  Gill  anemometer.  Figures  6  and  7 
show  that  the  sonic  anemometer  will  probably  under- 
estimate heat  and  momentum  fluxes  by  about  5%  to  30%, 
and  figure  8  shows  that  the  fluxes  associated  with  the 
three  chemical  species  will  be  underestimated  for  wind 
speeds  below  about  7  m/s  between  20%  or  30%  to  60%, 
depending  upon  stability.  At  wind  speeds  above  7  m/s, 
there  is  little  improvement  in  the  estimates  of  the  fluxes 
of  03,  N02,  and  NOx  with  the  sonic  anemometer  over 
the  Gill  anemometer;  this  is  largely  a  reflection  of  the 
nature  of  the  chemical  sensors  and  their  characteristics 
and  is  quite  independent  of  the  nature  of  the  vertical 
velocity  sensor. 

The  dotted  curves  of  figure  8  also  show  that  the  mis- 
matching of  these  two  types  of  sensors  can  introduce  an 
additional  5%  to  10%  underestimation  of  the  chemical 
fluxes,  depending  upon  atmospheric  stability  and  wind 
speed.  A  similar  computation  for  the  heat  flux  measure- 
ments showed  that  the  mismatching  errors  associated 
with  the  temperature  sensor  and  the  vertical  sonic 
anemometer  were  negligible. 

It  is  also  important  to  note  that  at  higher  wind  speeds 
the  mismatch  between  the  sonic  anemometer  and  the 
chemical  sensors  is  actually  worse  than  it  is  between  the 
Gill  propeller  anemometer  and  these  chemical  sensors. 
Consequently,  there  is  a  greater  loss  of  flux  at  higher 
wind  speeds  when  using  a  sonic  anemometer  to  measure 
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Figure  6.— Same  as  figure  3  except  for  the  use  of  a  sonic 
anemometer  rather  than  the  Gill  propeller  anemometer. 
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Figure  7.— Same  as  figure  4  except  for  the  use  of  the  3-axis  sonic 
anemometer  rather  than  the  3-axis  Gill  propeller  anemometer. 
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Figure  8.— Same  as  figure  5  except  for  the  use  of  the  sonic 
anemometer  rather  than  the  Gill  propeller  anemometer. 

the  fluctuating  component  of  the  vertical  velocity  than 
when  using  a  Gill  propeller  anemometer.  For  correcting 
our  flux  measurements,  we  assume  the  instruments  are 
mismatched  and  use  the  dashed  curves  in  figures  3,  5, 
and  8. 

Although  it  is  important  to  use  the  results  from  these 
figures  to  correct  turbulence  data  obtained  with  the  sonic 
anemometer,  it  is  also  important  to  correct  the  data  for 
shadowing  effects.  Shadowing  will  cause  the  measured 
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velocity  component  to  be  underestimated  because  of  par- 
tial shadowing  of  the  acoustic  path  by  the  tranducers. 
We  plan  to  use  the  algorithm  being  developed  by  Con- 
klin  et  al.  (1988)  for  this  correction  to  the  3-axis  sonic 
anemometer. 

In  addition  to  comparing  and  correcting  flux  estimates 
derived  with  the  Gill  and  sonic  anemometers,  other 
aspects  of  sensor  performance  were  considered.  One  was 
the  relative  importance  of  the  volume  averaging  effects 
that  the  intake  tube  of  the  chemical  sensors  can  have 
upon  the  estimates  of  the  flux  loss.  These  volume  averag- 
ing effects  arise  because  measurements  of  concentration 
are  not  made  at  a  point,  but  rather  such  measurements 
represent  the  average  concentration  of  a  finite  volume 
of  air  that  has  been  drawn  into  the  chemical  sensors 
through  the  Teflon  tube.  To  assess  these  effects,  esti- 
mates of  the  03,  N02,  and  NOx  flux  losses  with  and 
without  the  response  function  appropriate  to  volume 
averaging  (appendix  C)  were  needed  for  both  types  of 
anemometers.  Volume  averaging  effects  were  shown  to 
be  negligible  in  both  cases.  Therefore,  the  major  reason 
for  the  significant  losses  of  03,  N02,  and  NOx  fluxes 
(figs.  5  and  8)  is  largely  the  slow  response  times  of  these 
sensors. 


Variance  Loss  Corrections 

In  addition  to  the  flux  loss  correction,  the  variance  loss 
for  the  sensors  was  also  estimated.  Appendix  B  gives  the 
details  of  the  computations  used  in  the  following 
discussion. 

The  variance  loss  estimate  for  temperature,  vertical 
velocity  (for  both  the  Gill  and  sonic  anemometers)  and 
the  chemical  species  are  presented  in  figures  9-12.  The 
conclusions  that  can  be  drawn  from  these  figures  are  the 
same  as  with  the  flux  loss  estimates;  namely,  (1)  that  the 
sonic  anemometer  generally  represents  an  improvement 
over  the  Gill  propeller  anemometer,  and  (2)  that  the  slow 
response  of  the  chemical  sensors  can  cause  significant 
losses  to  the  variance  measurements,  especially  for  stable 
conditions.  However,  there  is  one  important  difference 
between  these  variance  loss  estimates  and  the  preceding 
flux  loss  estimates.  In  many  cases,  the  variance  loss  is 
explicitly  dependent  upon  the  stability  parameter  z/L,  as 
discussed  in  more  detail  in  appendix  B.  This  is  always 
true  whenever  the  atmosphere  is  stable,  and  it  is  true  for 
the  vertical  velocity  variance  for  unstable  conditions  as 
well.  For  the  flux  loss  estimate,  on  the  other  hand,  it  is 
sufficient  to  know  only  whether  the  atmosphere  is  stable 
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Figure  9.— Estimates  of  temperature  variance  loss  with  the  fast- 
response  platinum  wire  resistance  thermometer  as  a  function  of 
wind  speed  and  atmospheric  stability. 
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Figure  11.— Estimates  of  vertical  velocity  variance  loss  with  a 
sonic  anemometer  as  a  function  of  wind  speed  and  atmospheric 
stability. 
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Figure  10.— Estimates  of  vertical  velocity  variance  loss  with  the  Gill 
propeller  anemometer  as  a  function  of  wind  speed  and  at- 
mospheric stability. 
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Figure  12.— Estimates  of  variance  loss  from  03,  N02,  and  NOx 
measurements  as  a  function  of  wind  speed  and  atmospheric 
stability. 
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or  unstable  (equivalent  to  knowing  only  the  sign  of  z/L 
or  the  direction  of  the  heat  flux).  Therefore,  in  order  to 
correct  for  any  loss  of  variance,  we  must  first  correct 
for  any  loss  of  heat  and  momentum  flux  using  the  esti- 
mates derived  above;  next  we  compute  a  value  for  z/L 
based  on  the  heat  and  momentum  fluxes;  and  then  we 
make  an  estimate  of  variance  loss  for  that  value  of  z/L. 
Obviously,  estimating  the  variance  loss  is  considerably 
more  complicated  than  for  estimating  flux  losses. 

Finally,  an  important  caveat  must  be  stated  concern- 
ing all  estimates  of  the  flux  or  variance  losses.  There  is 
some  uncertainty  in  the  atmospheric  spectra  and  co- 
spectra  used  for  these  computations  [see  Moore  (1986) 
for  further  discussion].  Although  this  uncertainty  is  not 
serious  for  small  estimates  of  loss  (less  than  about  30%), 
it  can  become  more  important  at  higher  estimates  of  loss, 
which  can  make  the  correction  less  reliable  whenever 
the  flux  loss  is  large. 

Additional  Sources  of  Errors 

In  addition  to  the  errors  associated  with  instrument 
response,  size,  and  siting  geometry,  outlined  above,  it 
is  equally  important  to  be  aware  of  and  account  for  other 
important  measurement  limitations  and  problems  asso- 
ciated with  the  eddy  correlation  technique.  For  exam- 
ple, both  heat  and  moisture  fluxes  can  significantly 
influence  flux  measurements  (Webb  et  al.  1980),  as  can 
digitization  noise  (Kristensen  and  Kirkegaard  1987),  ran- 
dom uncorrected  noise  (Lenschow  and  Kristensen  1985), 
flow  distortion  effects  (Wyngaard  1988,  Hogstrom  1988), 
and  site  inadequacies  such  as  insufficient  fetch  and  lack 
of  surface  uniformity. 

Corrections  to  dry  deposition  fluxes  due  to  heat  and 
moisture  fluxes  arise  because  of  the  atmospheric  densi- 
ty variation  associated  primarily  with  the  heat  flux  and 
secondarily  with  the  moisture  flux.  For  our  site  the 
humidity  is  low  and  much  of  the  time  evapotranspira- 
tion  (moisture  flux)  is  not  large;  therefore,  we  will  use 
the  following  expression  for  correcting  the  measured 
fluxes  of  03  and  N02: 

F  =  Fm  +  (pc  wTW  [1] 

where  Fm  is  the  measured  flux  after  the  instrument  cor- 
rections outlined  in  the  previous  section  have  been  ap- 
plied, pc  is  the  mean  concentration  of  the  constituent 
(03  or  N02),  T  is  the  mean  ambient  temperature,  and 
wT'  is  the  heat  flux  also  measured  after  the  instru- 
ment corrections  have  been  applied.  Our  experience  has 
shown  that  the  correction  term  (pc  wT')/T  in  this  equa- 
tion is  typically  between  1%  and  20%  of  the  measured 
ozone  flux  and  is  frequently  significantly  higher  for  the 
measured  N02  fluxes.  The  NOx  sensoi,  on  the  other 
hand,  does  not  require  this  correction.  Fundamentally 
this  sensor  measures  the  mixing  ratio  rather  than  a  con- 
centration as  do  the  ozone  and  N02  sensors.  For  meas- 
urements of  fluxes  based  upon  mixing  ratios,  Webb  et 
al.  (1980)  have  shown  that  no  correction  need  be  applied. 


In  the  case  of  variance  measurements,  digitization 
noise  can  also  influence  measurements  (Kristensen  and 
Kirkegaard  1987).  Although  a  correction  for  this  effect 
was  included  in  our  data  analysis,  usually  it  was 
negligible. 

Uncorrelated  noise  (Lenschow  and  Kristensen  1985) 
can  also  impact  flux  measurements.  This  is  particularly 
true  whenever  the  instruments  are  event  counters.  Our 
chemical  sensors  measure  a  current  rather  than  count- 
ing photons,  and  the  concentrations  being  measured  are 
sufficiently  high  that  this  type  of  uncorrelated  noise  error 
does  not  contribute  significantly  to  our  errors. 

Flow  distortion  effects  caused  by  sensor  geometry  and 
their  positioning  on  towers  can  also  influence  observa- 
tions significantly  (Wyngaard  1988,  Hogstrom  1988).  Our 
measurements  are  undoubtedly  influenced  to  some 
extent  by  these  effects;  however,  drawing  air  samples 
into  the  shelter  with  an  intake  tube  minimizes  the 
amount  of  equipment  placed  on  the  towers,  which 
should  considerably  reduce  flow  distortion  effects. 

Finally,  patchiness  in  the  local  vegetation  cover  and 
inadequate  fetch  can  also  interfere  with  flux  measure- 
ments; but,  these  site  limitations  do  not  appear  at  pres- 
ent to  have  had  any  significant  impact  on  our  data.  (We 
say  this  because,  as  will  be  shown  in  the  next  section, 
our  data  agree  with  previously  and  independently  made 
observations.)  However,  eddy  correlation  measurements 
of  moisture  flux  and  independent  measurements  of  the 
surface  energy  balance  are  planned  and  should  provide 
more  insight  into  how  much  site  limitations  (and  flow 
distortion  as  well)  can  influence  the  data. 

DATA  RESULTS  AND  ANALYSIS 

Sample  results  from  our  Pawnee  program  were  used 
to  assess  the  utility  and  performances  of  the  different 
meteorological  and  chemical  sensors.  To  examine  nor- 
malized variances  of  temperature  and  vertical  velocity 
and  the  correlation  coefficient  between  temperature  and 
vertical  velocity,  the  Monin-Obukhov  similarity  theory 
was  used.  Similar  computations  were  performed  for  the 
ozone,  N02,  and  NOx  fluxes.  Deposition  velocities  and 
surface  resistances  derived  from  our  system  at  the 
Pawnee  site  were  compared  with  similar  measurements 
made  at  other  sites  for  these  three  acid  species. 

For  this  analysis,  we  selected  a  2-week  period  begin- 
ning on  June  30,  1987,  at  9:00  a.m.  MST,  which  made 
available  a  set  of  more  than  600  potentially  useful  data 
points  (fluxes,  variances,  etc.)  at  half-hour  intervals.  As 
an  example  of  these  data,  figures  13,  14,  and  15  give  the 
corrected  heat,  momentum,  and  ozone  fluxes,  respective- 
ly, for  the  week  beginning  with  June  30.  These  data  have 
not  been  smoothed  or  checked  for  bad  data  points 
(discussed  below),  and  breaks  in  the  curves  occur  when- 
ever data  were  not  recorded.  However,  not  all  these  data 
are  useful,  and  the  following  criteria  were  developed  for 
limiting  the  analysis  to  the  best  data. 

1.  Data  taken  when  the  wind  direction  was  between 
150°  and  225°  were  eliminated.  This  precluded  any  in- 
terference by  the  tower  and  the  instrument  shelter  on 
the  measurements. 
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2.  Data  taken  when  the  horizontal  wind  speed  was  less 
than  about  1.9  m/s  were  also  eliminated.  This  criterion 
was  necessary  to  ensure  that  the  anemometry  would 
resolve  the  momentum  flux  adequately.  The  exact  value 
of  this  cutoff  wind  speed  was  determined  by  trial  and 

0  5 
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Figure  13.— Corrected  heat  flux  for  week  beginning  June  30 
(day  181). 
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14.— Corrected  momentum  flux  for  week  beginning  June  30 
(day  181). 
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Figure  15.— Corrected  ozone  flux  for  week  beginning  June  30 
(day  181). 


error  and  resulted  from  an  analysis  of  the  roughness 
length  appropriate  to  our  site,  as  is  discussed  more  fully 
later. 

3.  All  momentum  flux  measurements  larger  than 
-0.0025  m2/s2  were  also  eliminated.  The  purpose  here 
was  to  remove  all  positive  momentum  fluxes  and  all 
small  negative  stresses  that  are  below  the  instrumenta- 
tion's ability  to  resolve.  Contrary  to  expectations,  the 
observed  momentum  flux  is  sometimes  positive,  possibly 
because  of  the  brief  and  intense  "bursts,"  which  in  the 
atmospheric  boundary  layer  have  been  observed  to 
reverse  the  sign  of  the  momentum  flux  (Haugen  et  al. 
1971).  These  positive  momentum  fluxes  were  removed 
because  all  computations  involving  Monin-Obukhov 
similarity  theory  and  surface  roughness  lengths  require 
negative  Reynolds  stress  only. 

4.  All  time  series  (temperature,  velocities,  chemical 
concentration,  fluxes,  etc.)  were  visually  checked  and  all 
obviously  noisy  data  points  were  removed.  These  data 
were  usually  in  the  form  of  individual  spikes  in  the  data 
and  were  simply  unrealistic  on  physical  grounds. 
(Examples  of  this  type  of  data  would  include  ambient 
temperature  measurements  that  exceed  the  boiling  point 
of  water.) 

5.  For  the  analysis  of  the  N02  and  NOx  data,  our 
criterion  required  that  the  observed  NOx  concentration 
exceed  the  observed  NOz  concentration  (otherwise  the 
data  from  both  sensors  were  eliminated).  In  a  few  in- 
stances, the  observed  N02  concentration  did  exceed  the 
observed  NOx  concentration.  In  some  cases,  this  was 
due  to  low  or  nonexistent  flow  rates  of  the  chemilumines- 
cent  reactant  through  the  reaction  chamber;  in  other 
cases,  however,  this  discrepancy  was  not  fully  under- 
stood but  may  have  been  due  to  problems  in  the  calibra- 
tion of  these  two  sensors. 

6.  All  data  referring  to  the  fast  response  chemical  sen- 
sors taken  during  periods  of  calibration  of  the  chemical 
instruments  were  eliminated.  During  a  calibration  period 
(usually  twice  weekly),  the  instruments  were  manually 
being  adjusted  and  tested;  hence  the  data  recorded  by 
the  computer  system  during  this  time  will  not  reflect  am- 
bient conditions.  A  calibration  of  these  instruments 
typically  requires  between  30  and  60  minutes. 

Employing  these  six  criteria  reduced  the  total  number 
of  acceptable  data  points  from  a  potential  of  more  than 
600  to  between  200  and  300,  depending  on  what  param- 
eter was  being  analyzed. 


Site  Roughness  and  Monin-Obukhov  Similarity  Theory 

Wind  speed  above  a  surface  is  influenced  by  the 
momentum  flux,  the  surface  roughness  length,  and  the 
atmospheric  stability.  The  relationship  between  these 
variables  is  defined  as  follows: 


u  =  ^  [ln(z/zQ)  -  tfm(z/L)] 


[2] 


where  \pm  is  the  atmospheric  stability  influence  function 
for  momentum,  z0  is  the  surface  roughness  length,  z  is 
the  instrument  height  (6  m),  k  is  the  von  Karman  con- 
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stant,  and  u,  is  the  friction  velocity  defined  as  (-u'w')1/2 
where  -u'w'  is  the  Reynolds  stress. 

The  influence  function,  \pm,  is  generally  expressed  as 
follows  (Paulson  1970,  Businger  et  al.  1971): 


2  lnp^)  +  lnp4r^)  -  tan^x  +  7r/2 
z/L<0  (unstable,  neutral) 

-5  z/L       z/L>0  (stable) 


[3] 


where  x  =  (1  -  16  z/L)1'4. 

Therefore,  using  our  measurements  of  heat  and 
momentum  fluxes,  wind  speed,  and  temperature  and 
employing  the  corrections  outlined  in  the  previous  sec- 
tion, we  could  compute  z0  with  the  above  relationships. 
Figure  16  shows  the  computed  values  of  roughness 
length  plotted  against  wind  speed  for  our  data  set  for  all 
wind  speeds  greater  than  1.9  m/s.  Except  for  one  outlier 
(z0  ~  45  cm,  u  ~3  m/s),  all  roughness  lengths  are  between 
about  0.5  and  30  cm.  The  reason  1.9  m/s  was  chosen  as 
the  lower  cutoff  speed  for  acceptable  measurements  is 
that  all  other  outliers  were  associated  with  these  lower 
wind  speeds  and  not  with  wind  direction.  Therefore,  it 
was  concluded  that  our  data  system  could  not  resolve 
momentum  flux  at  wind  speeds  below  about  1.9  m/s. 

The  average  z0  computed  from  the  data  shown  in 
figure  16  is  4.2  cm  ±  5.0  cm,  a  value  in  good  agreement 
with  the  4  cm  found  by  Korrell  et  al.  (1982)  for  most  wind 
directions  at  the  Boulder  (Colorado)  Meteorological 
Tower  which  is  located  about  100  km  southwest  of  our 
site  in  similar  terrain  with  similar  plant  cover.  However, 
our  results,  unlike  theirs,  did  not  show  any  correlation 
between  zD  and  wind  direction,  although  with  our  6-m 
measurement  height  our  z0  values  are  probably  influ- 
enced by  local  terrain  features. 

In  addition  to  the  large  variability  in  our  z0  values,  it 
is  also  important  to  examine  other  effects  that  may  lend 
further  uncertainty  to  defining  a  z0  value  for  our  site. 
Both  the  form  chosen  for  \j/m  and  any  systematic 
underestimation  of  the  momentum  flux  can  cause  zQ 
values  to  differ  significantly  from  those  shown  in  figure 
16.  For  example,  the  results  presented  by  Carl  et  al.  (1973) 
yield  a  different  form  for  ^m  appropriate  to  unstable  at- 
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Figure  16.— Roughness  length  vs.  wind  speed  for  the  2-week  period 
June  30-July  14. 


mospheric  conditions;  and,  Hogstrom  (1988)  found, 
when  correcting  [3]  for  flow  distortion  errors,  that  for 
unstable  condition  x  =  (1-15.2  z/L)1'4  and  for  stable  con- 
ditions ^m  (z/L)  =  -4.8  z/L.  These  two  corrections 
resulted  in  a  better  representation  of  ipm  (z/L)  than  [3]. 
The  results  of  Carl  et  al.  (1973)  produce  the  following 
expression  for  unstable  ^m: 

^m(z/L)  =  3/2  ln(1  +  x  +  x2)  _^3lan-f^±)  +  tt/vT 
3  3 


where  x  =  (1-15  z/L)1'3. 


z/L>0 


[4] 


Using  [4]  instead  of  [3]  yields  z0  values  that  are  about 
40%  lower  for  unstable  atmospheric  conditions  than 
those  shown  in  figure  16,  while  Hbgstrom's  formulation 
of  [3]  produced  essentially  negligible  changes  in  the  zQ 
values. 

Likewise,  systematic  underestimation  of  the  momen- 
tum flux  can  also  impact  our  z0  computations.  Such 
systematic  errors  can  arise  due  to  poor  vertical  align- 
ment, degradation  of  the  bearings  in  the  Gill  propeller, 
etc.  We  repeated  our  z0  computations  assuming  that  the 
momentum  stress  was  consistently  underestimated  by 
25%  and  found  that  all  z0  values  were  50%  higher  than 
those  shown  in  figure  16.  Given  these  results  and  those 
concerning  ^m  and  the  variability  in  zD  in  general,  we 
take  z0  =  4  cm  for  our  site,  but  the  uncertainty 
associated  with  this  value  is  large.  Finally,  it  should  be 
noted  that  z0  should  vary  with  time  of  year  and  the 
nature  of  the  plant  cover  or  snow  cover.  Therefore,  our 
z0  value  may  be  valid  only  for  the  present  data  set. 

Monin-Obukhov  similarity,  in  addition  to  providing 
some  insight  into  the  influence  function,  \pm  ,  for  our  z0 
correlations,  also  predicts  the  functional  relationships 
between  z/L  and  the  normalized  variances  of  vertical 
velocity,  aw2/u*2,  and  temperature,  (-z/L)<7T2u,2/(wT')2. 
Because  these  two  normalized  variances  are  combina- 
tions of  various  quantities  measured  by  our  system,  com- 
puting them  with  our  data  provides  (1)  a  method  to  assess 
the  quality  or  information  content  of  our  data  set  and 
(2)  a  means  for  checking  the  correction  factors  derived 
in  the  previous  section. 

Figure  17  is  a  plot  of  <jw2/u*2  vs.  z/L  as  computed  with 
our  data  and  as  expected  from  similarity  theory  [solid 


z/L 


Figure  17.— Vertical  velocity  variance  aw2  normalized  by  u  ,2  vs.  z/L. 
Observed  data  is  given  by  similarity  theory  prediction  is  given 
by  the  solid  line. 
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line  =  1.7(1-3  z/L)2'3,  from  Panofsky  et  al.  (1977)].  Al- 
though the  individual  data  points  obviously  show  signifi- 
cant variability,  the  data  in  some  average  sense  appear 
to  follow  the  expected  curve.  This  suggests  that  the  cor- 
rection factors  discussed  in  the  previous  section  are 
realistic  overall,  but  that  they  may  be  in  error  for  indi- 
vidual points.  This  may  arise  in  part  from  the  assump- 
tion that  the  response  length  of  the  vertical  Gill  propeller 
is  constant  (appendix  B);  whereas,  it  is  in  fact  a  func- 
tion of  the  wind  velocity's  angle  of  attack,  which  may 
vary  from  data  point  to  data  point.  Some  portion  of  the 
scatter  shown  in  figure  17  is  also  probably  due  to  errors 
in  the  measurement  of  the  momentum  flux.  This  is  a  par- 
ticularly demanding  measurement  which  is  quite  sen- 
sitive to  even  small  errors  in  alignment  and  other  effects. 
Finally,  during  stable  atmospheric  conditions  (z/L>0), 
these  difficulties  are  further  compounded  because  (1)  the 
fluxes  and  variances  will  usually  be  much  smaller  in 
magnitude  than  for  unstable  atmospheric  conditions, 
placing  greater  demands  on  the  instrumentation,  and  (2) 
the  correction  factors  become  less  certain  for  stable  con- 
ditions because  the  spectra  and  cospectra  are  themselves 
less  well  defined  in  stable  conditions.  

Figure  18  is  a  plot  of  (-z/L)aT2u*2/(wT')2  vs.  z/L  as 
computed  with  our  data  and  as  expected  from  similari- 
ty theory  (solid  line  =  0.92(-z/L]2'3,  from  Wyngaard  et 
al.  (1971).  The  agreement  between  our  observations  and 
similarity  theory  is  better  than  in  the  previous  figures 
and  is  actually  surprisingly  good  except  for  a  few  points 
between  z/L  =  0  and  z/L  =  -1.5.  The  reduction  in  scat- 
ter over  the  results  in  the  previous  figure  may  be  due 
to  the  use  of  a  relatively  high-performance  temperature 
sensor.  Certainly  fewer  assumptions  were  involved  in 
estimating  variance  and  flux  losses  with  this  sensor  than 
with  the  Gill  propeller  anemometer. 

Figure  19  is  a  plot  of  the  correlation  coefficient,  rwT 
between  temperature  and  vertical  velocity  vs.  z/L.  The 
correlation  coefficient  is  defined  as  wT7(ctw<jt).  This 
plot  shows  that  in  unstable  conditions  rwT~  +  0.55,  and 
for  stable  conditions  rwT— 0.15.  However,  the  data 
possess  a  large  amount  of  variability.  For  unstable  con- 
ditions, our  value  of  rwT  is  in  good  agreement  with 
Haugen  et  al.  (1971).  For  stable  conditions,  however, 
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Figure  19.— Correlation  coefficient  between  the  fluctuations  of  ver- 
tical velocity  and  temperature  vs.  z/L  for  the  2-week  period  June 
30-July  14. 


Haugen  et  al.  (1971)  suggest  that  rwT  -0.3.  Again,  this 
discrepancy  could  be  the  result  of  the  many  possible  er- 
rors associated  with  these  measurements  during  stable 
conditions. 


Ozone 

Up  to  now,  we  have  simply  assumed  that  ozone  is  trans- 
ported and  dispersed  by  atmospheric  turbulence  in  a 
manner  similar  to  temperature.  However,  comparing  the 
normalized  ozone  variance  (-z/L)a032  u*z/(w,o3')2  against 
z/L,  provides  an  important  consistency  check  of  this 
similarity  assumption.  Figure  20  is  the  same  as  figure 
18  except  for  a  different  range  of  z/L  values  and  the  use 
of  the  ozone  data  rather  than  the  temperature  data. 
Although  there  is  somewhat  more  scatter  and  fewer  data 
points  in  figure  20  than  in  figure  18,  the  results  are  still 
quite  similar.  There  is  nothing,  therefore,  in  either  of 
these  figures  to  invalidate  our  assumptions  of  similarity 
between  these  two  atmospheric  scalars. 

The  correlation  coefficient  between  ozone  and  the  ver- 
tical velocity,  rwc,3,  can  also  provide  some  information 
concerning  the  utility  of  the  variance  method  of  esti- 
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Figure  18.— Normalized  temperature  variance  vs.  z/L  for  the  2-week 
period  June  30-July  14.  Observed  data  is  given  by  •;  similarity 
theory  prediction  is  given  by  the  solid  line. 
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Figure  20.— Normalized  ozone  variance  vs.  z/L  for  the  2-week  period 
June  30-July  14.  Observed  data  is  given  by  •;  similarity  theory 
prediction  is  given  by  the  solid  line. 
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mating  the  ozone  flux  (Businger  1986).  From  figure  21, 
which  gives  a  plot  of  rwc,3  vs.  z/L,  rw03  can  probably 
be  approximated  by  a  single  value  (-0.55)  reasonably 
well.  Although  again  these  data  show  a  lot  of  scatter,  they 
are  in  good  agreement  with  the  rwT  for  unstable  condi- 
tions. However,  for  stable  conditions  all  values  between 
-0.1  and  -0.6  seem  to  be  about  equally  represented  for 
rwc,3.  Thus,  it  is  not  possible  to  confirm  or  reject  the 
usefulness  of  the  variance  method. 

In  addition  to  examining  the  ozone  data  using  Monin- 
Obukhov  similarity  theory,  we  can  also  use  this  data  to 
derive  deposition  velocities  and  surface  resistances.  The 
deposition  velocity,  vd,  is  defined  as  the  negative  turbu- 
lent flux  divided  by  the  mean  concentration;  the  surface 
resistance,  rs,  is  defined  as  follows: 


r.=  ^-k^'ln(zfaJ-^'L)]-^3 


[5] 


where  \pc  is  the  influence  function  for  ozone.  Assuming 
that  \pc  is  equal  to  the  influence  function  for  heat,  we 
have  (Paulson  1970):  . 


iMz/L) 


2ln( 


l  +  x< 


-5  z/L 


z/l<:o 

z/L>0 


[6] 


where  x  =  (1-16  z/L)1'4.  From  the  definition  of  surface 
resistance,  [5],  we  see  that  it  is  a  residual  and  equal  to 
the  total  resistance,  rt  =  l/vd,  less  the  aerodynamic 
resistance,  ([ln(z/z0)  -  \pc  (z/L)]/ku*),  and  the  boundary 
layer  resistance,  (2.63/ku*).  Further  discussion  of  these 
latter  two  resistance  terms  can  be  found  in  Hicks  et  al. 
(1985). 

Figures  22  and  23  give,  respectively,  the  deposition 
velocities  and  the  surface  resistances  for  ozone  as  a  func- 
tion of  the  time  of  day  for  our  data  set.  The  mean  deposi- 
tion velocity  for  ozone  is  0.43  cm/s  and  the  range  of  our 
observed  deposition  velocities  (about  0.1  cm/s  to  about 
1.2  cm/s)  are  in  good  agreement  with  all  other  observa- 
tions (e.g.,  Wesely  et  al.  1978;  Sehmel  1980;  Delaney  et 
al.  1986).  The  surface  resistance  for  ozone  basically 
should  reflect  the  behavior  of  the  plant  cover  during  the 
course  of  a  diurnal  cycle,  which  is  roughly  what  figure 
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Figure  21.— Correlation  coefficient  between  the  fluctuations  of 
vertical  velocity  and  ozone  concentration  vs.  z/L  for  the  2-week 
period  June  30- July  14. 
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Figure  22.— Ozone  deposition  velocity  as  a  function  of  time  of  day 
for  the  2-week  period  June  30-July  14. 
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Figure  23.— Surface  resistance  to  ozone  uptake  as  a  function  of 
time  of  day  for  the  2-week  period  June  30-July  14. 

23  shows.  Although  there  is  a  lot  of  scatter  in  these  data, 
generally  higher  resistances  are  seen  at  night  and  lower 
resistance  in  the  daylight  hours,  with  a  minimum 
resistance  between  2  and  about  4  s/cm  between  about 
900  hours  and  1,600  hours  and  a  maximum  resistance 
between  about  10  to  20  s/cm  at  other  times  during  the 
day.  This  diurnal  cycle  of  surface  resistance,  as  por- 
trayed in  the  figure,  is  in  good  agreement  with  the  diur- 
nal cycle  of  stomatal  resistance  for  Bouteloua  gracilis 
(Monson  et  al.  1986),  the  dominant  plant  in  the  area. 


N02  and  NOx 

The  N02  and  NOx  data  were  analyzed  in  a  manner 
identical  to  that  discussed  above  for  ozone  and,  in- 
asmuch as  the  results  and  conclusions  apply  equally  to 
both  N02  and  NOx,  only  the  NO,  sensor  will  be  dis- 
cussed. The  results  for  the  N02  sensor  are  summarized 
in  the  next  four  figures. 

Figures  24  and  25  give  the  measured  N02  concentra- 
tion (in  ppb/volume)  and  corrected  N02  fluxes  for  the 
first  week  of  our  data  set.  Both  the  NO,  concentrations 
and  the  N02  fluxes  are  highly  variable.  The  positive 
N02  fluxes  in  figure  25  indicate  that  the  surface  itself 
may  be  a  source  for  N02,  as  has  in  fact  been  observed 
(Slemr  and  Seiler  1984,  Williams  et  al.  1987).  Further- 
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Figure  24.  — N02  concentration  as  a  function  of  time  during  the 
1-week  period  June  30-July  7. 
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Figure  25.— N02  flux  as  a  function  of  time  during  the  1-week  period 
June  30-July  7. 

more,  the  variability  in  the  N02  source  strength  at  the 
surface  has  also  been  observed  to  be  extremely  variable 
in  both  time  and  space.  For  example,  Slemr  and  Seiler 
(1984)  reported  variation  in  source  strength  of  over 
100-fold.  To  further  complicate  the  interpretation  of  our 
data,  natural  NOx  emitted  from  soils  is  oxidized  in  the 
presence  of  ozone  and  ultraviolet  radiation  to  N02. 
When  this  occurs  (which  it  certainly  must  at  our  site), 
Fitzjarrald  and  Lenschow  (1983)  and  Lenschow  and 
Delany  (1987)  have  shown  that  the  measured  fluxes  can 
be  in  error  by  significant  amounts  (about  30%  at  our 
measuring  height  of  6  m). 

Figure  26  gives  rWN02  as  a  function  of  z/L  and  figure 
27  gives  the  N02  deposition  velocity  for  N02  as  a  func- 
tion of  the  time  of  day.  In  the  first  of  these  two  figures, 
there  is  no  clear  evidence  of  any  correlation  between  the 
vertical  velocity  fluctuations  and  the  N02  concentration 
fluctuations.  The  deposition  velocities  clearly  reflect  the 
variability  in  the  N02  flux  itself. 

Although  our  results  are  inconclusive  for  N02,  it 
should  be  noted  that  our  observations  are  very  similar 
to  other  observations  of  NOz  and  NOx  measurements 
(e.g.,  Delaney  et  al.  1986,  Hicks  et  al.  1986).  In  general, 
until  more  is  understood  quantitatively  about  the  role  of 
generation  and  destruction  of  NOx  and  N02  at  the  sur- 
face and  within  the  soil,  it  will  be  difficult  to  make  fur- 
ther progress  on  interpreting  similar  measurements. 
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Figure  26.— Correlation  coefficient  between  the  fluctuations  of 
vertical  velocity  and  N02  concentration  vs.  z/L. 
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Figure  27.  — N02  deposition  velocity  as  a  function  of  time  of  day 
for  the  2-week  period  June  30-July  14. 

Furthermore,  it  will  probably  require  some  sort  of  a  soil 
nitrogen  cycle  model  coupled  with  a  chemical  and  at- 
mospheric transport  model  to  fully  interpret  our  results. 

SUMMARY 

Our  efforts  at  measuring  the  dry  deposition  fluxes  of 
ozone,  N02,  and  NOx  employing  the  eddy  correlation 
method  have  been  outlined.  Although  this  method  is 
demanding  and  complex,  at  present  it  offers  the  only 
direct  approach  available  for  making  these  measure- 
ments and,  therefore,  we  chose  this  technique  over  any 
other.  Although  an  attempt  was  made  to  account  for  as 
many  errors  associated  with  this  approach  as  possible, 
sources  of  errors  undoubtedly  remain.  Nevertheless,  the 
results  of  our  initial  analysis  suggest  that  our  eddy  cor- 
relation system  produces  realistic  and  valid  data  that 
agree  with  many  previous  and  independent  meteor- 
ological and  biological  measurements. 

Finally,  it  should  be  emphasized  that  this  report  is  part 
of  a  continuing  effort  at  the  Pawnee  Grasslands.  More 
detailed  analysis  of  the  data  gathered  to  date  is  planned, 
and  the  system  will  be  expanded  to  include  moisture  flux 
measurements  using  eddy  correlation  methods.  Further- 
more, long-term  comparisons  of  the  eddy  correlation 
fluxes  of  heat  and  moisture  with  other  independently 
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measured  terms  of  the  surface  energy  balance  are  also 
planned.  Hopefully,  these  and  other  future  plans  should 
further  demonstrate  the  utility  and  validity  of  our  system 
and  the  quality  of  our  site. 
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APPENDIX  A.  Description  of  Fast  Response  Chemical 
Sensors1 

An  earlier  version  of  this  research  project  was  de- 
signed to  measure  the  dry  deposition  of  N02,  NOx 
(principally  NO  +  N02  +  peroxyacetyl  nitrate  [PAN]), 
ozone,  and  SOx  to  grassland  during  winter  and  summer 
months.  However,  the  SOx  sensor  was  not  available  dur- 
ing the  field  season  because  it  was  still  undergoing  fur- 
ther development  in  the  laboratory  to  improve  sensitivity 
and  response  time.  This  chemical  analysis  project  in- 
cluded the  installation,  operation,  and  maintenance  of 
the  intake  system,  detectors,  and  the  zero  and  calibra- 
tion system. 

The  instruments  were  housed  in  the  shelter,  and  air 
was  brought  from  the  tower  to  the  shelter  using  a  high- 
volume  sample  pump.  The  arrangement  avoided  in- 
troducing extraneous  turbulence  to  the  mean  wind  flow 
and  kept  the  instruments  in  a  temperature-controlled  en- 
vironment. The  intake  system  consisted  of  a  26.5-m-long, 
5/8-inch-i.d.  Teflon  tube  through  which  air  was  drawn 
at  ca.  160  1/min  using  a  high-volume  air  pump  (Fuji  ring 
compressor).  The  flow  rate  decreased  with  time  (table 
Al).  The  chemical  sensors  sampled  from  the  main  air 
stream  and  returned  their  exhaust  downstream  of  the 
sampling  manifold  to  keep  the  pressure  drop  across  the 
instruments  to  a  minimum.  The  interference  caused  by 
adsorption/desorption  was  minimized  by  using  large- 
bore  Teflon  tubing.  Although  the  length  of  sample  tube 
slowed  the  response  of  the  chemical  instruments  and  cut 
off  some  of  the  high-frequency  flux  data,  the  need  for 
a  complex  temperature  correction  was  eliminated  by 
housing  the  instruments  in  a  temperature-controlled 
shelter.  Table  A2  presents  a  summary  of  the  chemical 
instruments  and  their  specifications. 


Chemiluminescent  Ambient  Air  Monitors  (CAAM) 

The  same  instrument  design  was  adapted  to  measure 
nitrogen  dioxide  (Wendel  et  al.  1983)  and  ozone  (Ray  et 
al.  1986)  by  changing  the  fluorescent  reagent  used.  The 
sample  flowed  through  a  detection  cell  and  over  a  filter 
pad  wetted  with  the  fluorescent  reagent  (see  fig.  Al).  A 
blue-sensitive  photomultiplier  tube  (PMT)  monitored  the 
light  produced.  The  fluid  and  air  exited  at  the  bottom 
of  the  reaction  cell  and  were  separated  by  gravity  in  the 
catch  reservoir.  The  fluid  flow  rate  was  1  ml/min  and 
the  air  flow  rate  was  10  1/min.  The  response  of  these  in- 
struments was  pressure-dependent. 

For  the  detection  of  nitrogen  dioxide  a  basic  solution 
of  luminol  (5-amino-2,3-dihydro-l,4-phthalazinedione) 
was  used  (Maeda  et  al.  1980).  The  solution  was  com- 
prised of  0.05  M  NaOH,  0.1  M  Na2S03,  2  x  10  4  M 
luminol,  and  0.05%  vol/vol  t-butanol.  The  solution  was 


^Taken  from  the  Final  Report  to  the  USDA  Forest  Service, 
"Measurement  of  Dry  Deposition  of  SOx,  NO„,  and  N02  and  03," 
September  1987,  under  USFS  Contract  No.  28-K6-384.  Report  was 
written  by  D.W.  Stocker,  D.H.  Stedman,  B.A.  Watkins,  and  J. A.  Arm- 
strong of  the  University  of  Denver,  Chemistry  Department,  Denver, 
CO  80208. 
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Figure  A1.— Schematic  diagram  of  the  Chemiluminescent  Ambient 
Air  Monitor,  with  an  enlargement  of  the  detection  cell,  which  is 
machined  from  black  Delrin;  a  glass  window  seals  the  cell. 

used  once  and  then  discarded.  Operating  under  these 
conditions,  the  sensitivity  was  ca.  0.2  ^A/ppb  (parts  per 
billion  volumetric  mixing  ratio)  and  the  response  was 
linear  from  0.1  ppb  (fig.  A2).  The  response  of  the  instru- 
ment did  decay  with  time  after  a  fresh  solution  was  in- 
stalled and  the  filter  pad  in  the  detection  cell  changed, 
probably  because  of  degradation  of  the  filter  pad  and  ag- 
ing of  the  solution.  The  solution  was  replaced  approx- 
imately weekly  and  the  instrument  recalibrated  (table  Al 
gives  details  of  the  calibration  changes).  When  operated 
in  the  laboratory  the  instrument  responded  to  a  square 
wave  pulse  of  N02  introduced  to  the  intake  of  the  in- 
strument within  0.2  seconds.  Under  field  conditions,  the 
response  to  a  square  wave  pulse  of  clean  air  introduced 
to  the  system  at  the  tower  was  slower  (0.7  s)  because 
broadening  of  the  sample  occurred  in  the  intake  system 
(fig.  A3).  The  detection  limit  was  0.14  ppb. 

For  the  analysis  of  ozone  the  basic  luminol  solution 
was  exchanged  for  a  solution  of  eosine-Y  in  ethylene 
glycol.  The  dye  solution  was  recycled  after  each  batch 
had  passed  through  the  instrument.  The  ozone  detector 
was  calibrated  by  comparing  the  output  signal  with  that 
from  a  slow-response  UV  photometric  ozone  monitor 
(TECO  model  49A).  The  sensitivity  was  0.12  /xA/ppb  and 
the  detection  limit  was  <1  ppb  in  the  field.  In  the 
laboratory  the  response  of  the  instrument  was  7  Hz, 
which  decreased  under  field  conditions  to  1.4  Hz  be- 
cause of  the  broadening  that  occurred  when  sampling 
through  the  long  intake  tube  (fig.  A4). 
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Table  A1.— Dates  of  calibration  and  field  notes  for  the  chemical  sensors  operated  at  the  Pawnee 
Grasslands  field  site  (Dec.  1, 1986,  through  Feb.  2, 1987,  and  May  1, 1987,  through  July  27, 1987). 
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aMeasured  using  a  dry  gas  meter  attached  to  the  exhaust  of  the  large  sample  pump,  correct- 
ing for  the  flow  through  the  NOx  analyzer  that  is  not  returned  to  the  sample  inlet. 

bWherever  possible,  the  times  when  the  dye  solution  was  recycled  is  reported. 

cThe  calibration  factor  is  in  units  of  ppb/V,  the  calibrated  sensitivity  in  parts  per  billion  mix- 
ing ratio  per  volt  observed  by  the  data  acquisition  system  when  the  instrument  is  working  with 
a  full-scale  deflection  of  10  \iA. 

d10  nA  full-scale  deflection. 

eFresh  solution. 

'Converter  heater  repaired. 

^Solution  changed,  but  neither  old  nor  new  solutions  calibrated. 
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"Changed  target  pad. 

'Changed  full-scale  deflection  to  100  ^A. 

'Changed  zero  air  traps. 

kAir  pump  replaced. 

'Response  of  N0X  analyzer  has  changed  by  a  factor  of  100,  which  implies  an  electronic  rather 
than  a  chemical  fault  is  to  blame.  The  response  of  the  N0X  analyzer  is  also  very  noisy. 
mLuminol  flow  stopped. 
"No  calibration  of  old  solution. 
°No  luminol  left. 
POxygen  tank  empty. 

<JHigh  calibration,  signal  may  not  have  stabilized. 
'TECO  photometric  ozone  monitor  recalibrated. 
sNo  calibration  possible  as  N0X  analyzer  off  line. 
'Plugged  frit  in  ozone  line. 

"Repair  of  ozone  generator  in  Denver  University  laboratories. 
vChange  air  filters  in  main  sample  line,  system  is  off-line. 
^Reservoir  dry  before  refill. 
"Plugged  sample  intake. 
YBroken  air  pump. 

zZero  air  system  self  destructs,  trap  explodes. 
aaReplace  zero  air  trap. 

Table  A2.— Summary  of  the  characteristics  of  the  chemical  sensors  used  during  the  field  study. 


Chemical 
measured 


Instrument 
name 


Sensitivity 


Comments 


Total  N0X   NO,,  box     0.3  ^A/ppb 


NO, 


Ozone 
Ozone 


CAAM  2      0.2  HA/ppb 


Teco  49A 
CAAM  1 


0.2  mv/ppb 
0.12  fjA/ppb 


Requires  high-pressure  02  cylinders  for  03  production; 
large  vacuum  pumps;  cooled  red-sensitive  PMT;  high 
power  (240  v).  Fast  response  <1  sec  NOx  by  catalytic 
reduction  using  Mo  turnings  at  400°C.  Stable  calibration. 

Uses  a  basic  luminol  solution  which  cannot  be  recycled. 
Each  batch  of  reagent  requires  calibration.  Sensitivity  of 
solution  decreases  with  time.  Calibrate  and  check  target 
pad  weekly.  Response  time  <1  sec.  Portable  low  power 
required. 

EPA  standard.  Slow  response. 

Uses  eosine-Y  in  ethylene  glycol  which  can  be  recycled 
but  some  sensitivity  is  lost  with  use.  Response  time  <1 
sec.  Calibrate  using  UV  photometric  instrument.  Period 
of  unattended  use  limited  by  reservoir  size.  Target  pad 
checked  weekly.  Portable  low  power  required. 
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Figure  A2.— The  response  of  the  total  nitrogen  oxides  detector  and  the  nitrogen  dioxide  de- 
tector to  low  concentrations  of  nitrogen  dioxide. 
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Figure  A3.— Frequency  response  of  the  nitrogen  dioxide  detector: 
(a)  to  a  square  wave  pulse  of  nitrogen  dioxide  in  the  laboratory 
and  (b)  to  a  square  wave  pulse  of  zero  air  introduced  to  the  tower 
inlet  in  the  field. 

Total  NOx  Analyzer 

NO  reacts  with  ozone  forming  electronically  excited 
nitrogen  dioxide.  At  reduced  pressure  the  excited  nitro- 
gen dioxide  decays  to  the  ground  state,  emitting  a 
photon: 


NO  +  03 
NO^ 


>no: 


O. 


 >N02  +  hp  (640  nm) 

The  intensity  of  light  produced  is  proportional  to  the  con- 
centration of  NO,  forming  the  basis  for  the  total  NOx 
analyzer  (Fontijn  1970). 

The  higher  oxides  of  nitrogen  were  reduced  to  NO  by 
passing  the  sample  gas  through  a  hot  molybdenum  con- 
verter at  400°C.  The  principle  species  reduced  were 
NOz  and  PAN  (Dickerson  et  al.  1984).  Nitric  acid  is 
known  to  pass  through  the  length  of  the  sample  tubes 
to  reach  the  converter  housed  in  the  shelter,  but  it  has 
been  shown  in  the  past  to  have  a  response  time  larger 
than  120  s.  This  slow  response  is  presumed  to  arise  from 
a  "chromatographic"  effect  of  HN03  sticking  and 
releasing  from  the  intake  tubing.  The  response  time  for 
nitric  acid  is  too  slow  to  allow  correlation  with  changes 
in  vertical  wind  velocities  to  determine  HN03  fluxes. 
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Figure  A4.— Frequency  response  of  the  ozone  detector:  (a)  to  a 
square  wave  pulse  of  ozone  in  the  laboratory  and  (b)  to  a  square 
wave  pulse  of  zero  air  introduced  to  the  tower  inlet  in  the  field. 

The  reduced  sample  gas  was  mixed  with  excess  ozone 
prior  to  entering  the  reaction  chamber  in  front  of  a 
cooled  red-sensitive  photomultiplier  tube  (EMMI  9658). 
The  walls  of  the  reaction  chamber  were  gold-plated  to 
enhance  reflectivity  and  provide  an  inert  surface.  Ozone 
was  produced  by  passing  250  ml/min  of  02  through 
three  parallel  discharge  tubes;  a  water  bubbler  in  the 
ozone  line  decreases  the  background  signal  (Drummond 
and  Voltz  1982).  The  analyzer  is  shown  schematically 
in  figure  A5.  The  photomultiplier  tube  output  was  fed 
to  a  programmable  electrometer  (Keithley  18000-20). 

The  instrument  was  operated  with  a  chamber  vacuum 
of  5-10  torr  and  a  sample  flow  rate  of  approximately  4 
1/min.  Under  these  conditions  the  instrument  response 
was  ca.  0.3  nA/ppb  with  background  noise  equivalent  to 
a  detection  limit  of  ca.  0.07-0.2  ppb;  however,  at  times 
the  noise  was  greater  than  this.  The  sensitivity  of  the  in- 
strument was  constant  throughout  each  experimental 
period,  as  shown  in  table  Al.  The  change  in  sensitivity 
between  the  two  experimental  periods  was  due  to  a  fault 
in  the  electrometer.  The  instrument  was  less  sensitive 
when  ambient  air  was  used  for  ozone  generation,  which 
was  necessary  when  the  oxygen  cylinder  had  emptied. 
The  response  of  the  analyzer  to  NO,  was  linear  until 
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the  signal  to  noise  ratio  approached  unit  (fig.  A2).  When 
operated  in  the  laboratory,  the  instrument  responded  to 
a  square  wave  pulse  of  10  ppb  of  NO  within  0.3  s  (fig. 
A6).  However,  when  operated  in  the  field  and  sampling 
through  the  long  intake  tube,  the  instrument  responsed 
somewhat  more  slowly  at  0.9  s. 

Experiments  carried  out  by  heating  the  intake  gases 
to  the  analyzer  while  keeping  the  concentration  of  span 
gas  constant  showed  no  deviation  in  response;  therefore, 
the  response  of  the  instrument  was  not  pressure- 
dependent.  This  is  consistent  with  the  conclusions  of 
Steffenson  and  Stedman  (1974)  provided  constant- 
volume  vacuum  pumps  are  used. 

Zero  and  Calibration  Systems 

The  instruments  used  in  this  study  are  highly  sensitive 
and  consequently  are  susceptible  to  some  baseline  drift; 
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Figure  A5.— Schematic  diagram  of  the  total  oxides  of  nitrogen 
detector  with  a  sectional  view  of  the  reaction  chamber. 
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Figure  A6.— Frequency  response  of  the  total  NOx  analyzer:  (a)  to 
a  square  wave  pulse  of  10  ppb  NO  in  the  laboratory  and  (b)  to  a 
square  wave  pulse  of  zero  air  introduced  to  the  tower  inlet  in  the 
field. 

thus,  it  is  necessary  to  keep  track  of  the  response  of  these 
instruments  to  clean,  unpolluted  air. 

"Zero  air"  was  prepared  by  scrubbing  ambient  air  us- 
ing a  series  of  traps.  The  first  contained  ferrous  sulphate 
which  reduced  N02  to  NO  which  was  then  removed  by 
a  second  trap  filled  with  Purafil.  A  third  trap  filled  with 
activated  charcoal  removed  ozone.  A  compressor  was 
used  to  fill  a  30-gallon  tank  to  a  pressure  of  100  psi  with 
zero  air.  At  the  end  of  each  sampling  period  zero  air  was 
introduced  at  the  inlet  to  the  main  sample  line  through 
3/8-inch-i.d.  black  Teflon  tubing.  The  flow  of  air  to  the 
tower  intake,  ca.  170  1/m,  was  sufficient  to  overvent  the 
main  sample  flow  at  a  constant  flow  rate  for  approx- 
imately 3  minutes. 

The  system  was  calibrated  on  a  weekly  basis.  This  was 
achieved  by  adding  a  known  flow  of  span  gas  to  the  zero 
air  supply  while  the  zero  air  flow  rate  was  still  constant. 
The  NOx  analyzer  was  calibrated  using  a  span  gas  mix- 
ture of  10  ppm  NO  in  N2.  The  N02  CAAM  instrument 
was  calibrated  using  a  1000  ppm  N02  in  N2  standard; 
the  concentration  of  N02  reaching  the  instruments  was 
checked  against  the  NOx  analyzer. 

Future  Work 

As  soon  as  the  SOx  analyzer  is  functioning  in  a  reliable 
and  predictable  mode,  field  measurements  of  the  flux  of 
SOx  to  prairie  grassland  and  snow  should  be  perform- 
ed. Other  chemical  measurements  of  interest  include  the 
measurement  of  the  flux  of  NOy  (NOx  +  HN03)  using 
a  new  intake  system  (currently  under  development)  for 
the  chemiluminescent  NOx  analyzer  used  in  this  study. 
The  NOx  analyzer  in  its  present  configuration  is  capable 
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of  measuring  the  flux  of  NO,  this  along  with  the  concen- 
tration of  PAN  should  be  measured  to  determine  a 
nitrogen  flux  budget.  The  concentration  of  PAN  should 
be  determined  using  a  chromatographic  technique 
developed  in  these  laboratories  based  on  a  luminol  detec- 
tor (Burkhardt  et  al.  1987). 
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APPENDIX  B.  Eddy  Correlation  System  Response 
Characteristics 

This  appendix  outlines  the  method  we  used  to  correct 
the  flux  measurements  for  sensor  inadequacies  such  as 
slow  response  characteristics,  spatial  averaging  effects, 
and  sensor  mismatching.  The  method  is  similar  to 
Moore's  (1986)  in  that  it  uses  the  frequency-dependent 
instrument  response  functions  and  the  spectral  and 
cospectral  distribution  functions  associated  with  at- 
mospheric turbulent  variances  and  fluxes  to  derive  the 
correction  factors.  It  differs  from  Moore's  (1986)  method 
in  that  it  employs  chemical  sensors,  treats  aliasing  dif- 
ferently, and  includes  the  effects  of  mismatched  sensors. 

Specifically,  the  correction,  AF,  of  the  true  flux,  F,  of 
a  quantity  with  specific  density,  q,  is  estimated  as 
follows: 

AF  _  _  /rHwq(n)Cwq(n)dn  rB1] 
F  /rcwq(n)dn 

where  the  true  flux  F  =  J0  Cwq(n)dn,  AF/F  represents 
the  percentage  flux  loss,  and  Hwq(n)  is  the  total  system 
cospectral  response  function  associated  with  the  sensors 
for  the  vertical  wind  velocity,  w,  and  the  quantity,  q. 
Hwq(n)  is  the  product  of  the  separate  response  functions 
associated  with  the  data  logging  system,  sensor  frequen- 
cy response,  sensor  size,  sensor  separation,  and  sensor 
mismatch,  and  Cwq(n)  is  the  atmospheric  cospectrum  of 
w  and  q  as  a  function  of  frequency,  n. 

The  response  functions  appropriate  to  each  of  the 
various  instruments  and  the  data  logging  system  are 
described. 


Fast  Response  Chemical  Sensors  (03,  N02,  NOx) 

At  present  the  fast  response  chemical  sensors  are 
housed  in  an  air-conditioned  shelter  about  26.5  m  south 
of  the  meteorological  tower.  The  sampled  air  is  brought 
rapidly  into  the  shelter  through  a  long  Teflon  tube,  whose 
intake  is  positioned  on  the  tower  near  the  anemometers 
and  temperature  sensor.  The  response  times  of  the  sen- 
sors were  determined  in  situ  by  a  series  of  concentra- 
tion step  changes  that  involved  alternately  switching 
between  passing  ambient  air  and  passing  clean  air  (zero 


concentration)  through  the  intake  tube  and  into  the 
sensors. 

Figures  B1-B3  give  examples  of  the  data  gathered  dur- 
ing one  of  the  response  time  tests  and  table  Bl  gives  the 
results  of  the  tests.  The  tests  were  performed  on  three 
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Figure  B1.— Response  time  tests  for  fast  response  ozone  sensor, 
June  24,  1987. 
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Figure  B2.— Response  time  tests  for  fast  response  N02  sensor, 
July  24,  1987. 
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Figure  B3.— Response  time  tests  for  fast  response  NOx  sensor, 
July  24,  1987. 
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separate  occasions  and  usually  consisted  of  three  or  four 
positive  and  negative  step  changes.  For  a  negative  step 
change  (decreasing  from  ambient  to  zero  concentration), 
the  data  were  fit  with  a  decreasing  exponential,  e_t'T_, 
and  the  best  response  time,  r-,  was  found  using  a  least 
squares  minimization  technique.  (Here  t  represents  time.) 
For  a  positive  step  change  (increasing  from  zero  to 
ambient  concentration),  the  same  method  was  employed 
except  that  the  exponential,  e~t/r~,  was  replaced  by 
(1  -  e~t/T  +  ).  Finally,  only  data  between  20%  and  80%  of 
maximum  were  used  to  fit  the  curves.  The  ozone  sensor 
(fig.  Bl),  which  usually  provides  the  least  noisy  signal, 
has  a  response  time  of  about  0.7  s.  Although  its  response 
time  was  independent  of  the  magnitude  of  the  step 
change,  the  ozone  sensor  does  display  some  slight  asym- 
metry in  response  to  positive  and  negative  step  changes 
with  the  response  time  r-  being  slightly  shorter  than  t+  . 
This  asymmetry  is  also  observed  in  the  laboratory  and 
is  thought  to  be  related  to  adsorption  and  desorption  of 
ozone  by  the  intake  tube  walls.  The  response  time  of  0.7  s 
is  also  slower  (by  a  factor  of  about  2)  than  that  observed 
in  the  laboratory  (Ray  et  al.  1986  and  appendix  A).  Again, 
this  is  consistent  with  the  adsorption  and  desorption 
processes,  which  will  tend  to  lengthen  the  response  time 
by  adding  ozone  to  the  air  stream  as  the  concentration 
decreases  to  zero  and  by  removing  ozone  from  the  air 
stream  as  the  concentration  increases  to  ambient.  In  ad- 
dition, the  mixing  of  gas  samples  and  the  resulting 
smearing  of  the  high  frequency  changes  in  gas  concen- 
trations by  the  highly  turbulent  flows  inside  the  Teflon 

Table  B1.— Optimal  response  times  of  fast  response  chemical 
sensors.3 


Sensor  response  time  (seconds) 


Test  date 

o3 

(CAAM1) 

N02 
(CAAM2) 

NOx 

January  28,  1987 

T" 

0.65 

T.r 

0.66 

T 

0.60 

T+ 

1.02 

T+ 

0.94 

T+ 

1.07 

June  3,  1987 

T" 

0.62 

0.72 

T" 

0.63 

0.63 

T~ 

0.64 

0.48 

T  + 

0.79 

T  + 

0.83 

T  + 

0.71 

June  24,  1987 

T" 

0.58 

T" 

0.66 

T" 

0.67 

T" 

0.56 

0.95 

T+ 

0.76 

0.79 

0.93 

T+ 

0.87 

0.63 

0.84 

T+ 

0.83 

0.73 

0.84 

T+ 

0.84 

0.65 

0.79 

ar_  equals  response  time  to  step  change  from  ambient  to  clear 
air  and  t+  equals  response  time  to  step  change  from  clean  air  to 
ambient. 


tube  (Reynolds  number  about  104)  may  also  contribute 
to  this  lengthening  of  the  response  time  of  the  chemical 
sensors  over  what  has  been  observed  in  the  laboratory. 

The  response  times  of  the  N02  and  NOx  sensors  (figs. 
B2  and  B3,  respectively)  were  more  difficult  to  obtain 
because  the  sensors  are  noisier  than  the  ozone  sensor. 
For  the  N02  sensor  the  response  time  is  0.7  s  and  it  is 
0.9  s  for  the  NOx  sensor.  These  response  times  were 
found  by  the  same  method  as  the  ozone  sensor  except 
that  there  were  fewer  useful  data  sets  with  which  to  do 
the  curve  fitting.  In  the  case  of  the  NOx  sensor,  a  low 
pass  RC  analog  filter  with  a  0.04-s  time  constant  had  to 
be  incorporated  to  remove  the  high  frequency  noise  in 
the  signal.  These  two  sensors,  as  does  the  ozone  sensor, 
also  display  a  slight  asymmetry  in  their  response  to 
positive  and  negative  step  changes,  and  a  longer 
response  time  than  observed  in  the  laboratory  (Wendel 
et  al.  1983,  Wesely  et  al.  1982).  Again,  as  with  ozone,  ad- 
sorption and  desorption  of  these  gases  by  the  tube  walls 
and  the  turbulent  flows  inside  the  tube  are  thought  to 
be  largely  responsible  for  these  differences.  These  last 
two  figures  also  show  another  problem  peculiar  to  these 
two  sensors.  In  both  cases,  the  N02  and  NOx  concen- 
trations increase  above  ambient  just  prior  to  the  region 
where  the  concentrations  begin  to  decrease  to  zero.  The 
cause  for  this  pulse  of  gas  is  unknown,  but  it  may  be 
related  to  the  valves  that  switch  the  clean  air  flow  and 
the  ambient  air  flow.  For  the  response  time  tests  this 
problem  was  not  a  consideration. 

It  is  important  to  note  that,  in  addition  to  the  response 
times,  these  response  time  tests  also  give  the  lag  or  delay 
time  required  to  move  the  air  from  the  mouth  of  the  in- 
take tube  to  the  sensors.  The  delay  time  is  represented 
on  figures  B1-B3  as  the  time  from  input  of  clean  air  at 
the  mouth  of  the  intake  tube  (t~0)  to  the  time  at  which 
the  concentration  begins  to  decrease.  These  results  show 
that  the  lag  time  was  2.0  s  at  a  flow  rate  of  155  1/min. 
This  direct  method  of  determining  the  delay  time  was 
also  independently  verified  from  measurements  of  the 
flow  rate  and  the  length  and  cross-sectional  area  of  the 
intake  tube  (inside  diameter  of  intake  tube  =  5/8  inch). 
This  independent  estimate  of  the  delay  time  was  2.1s. 
For  computing  eddy  correlation  fluxes,  the  analog  sig- 
nals from  these  chemical  sensors  are  shifted  digitally  in 
time  by  2  s. 

Finally,  these  response  time  tests  also  show  that  the 
zero  concentration  of  these  gases  do  not  necessarily  cor- 
respond to  a  zero  voltage  reading  of  the  instruments. 
This  bias  in  the  data  signal  was  intentional  and  was  in- 
cluded in  data  calibration  and  reduction.  Further  details 
on  the  calibration  of  these  chemical  sensors  is  discussed 
in  appendix  A. 

With  these  response  times,  the  response  function  that 
describes  the  dynamic  frequency  response  of  these 
chemical  sensors  is  given  as  follows: 

Td(n)  =  l/[l  +  (27rnT)2]1'2  [B2] 

Here  we  assume  that  these  sensors  are  "first  order 
systems"  (Brock  and  Nicolaidis  1985). 
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In  addition  to  the  dynamic  frequency  response, 
another  important  question  that  needs  to  be  addressed 
is  that  of  spatial  averaging  effects  caused  by  sampling 
a  volume  of  air  as  it  is  drawn  through  the  intake  tube 
into  the  reaction  chambers.  This  effect,  which  is  likely 
to  be  important  whenever  the  sampled  volume  of  air  is 
large  and  the  sensor  response  is  sufficiently  fast,  results 
from  the  fact  that  a  measurement  of  concentration  is  not 
made  at  a  point,  but  rather  the  measurements  of  concen- 
tration represent  the  average  concentration  within  a 
finite  volume  of  air  that  has  been  drawn  into  the  sensor. 
These  spatial  averaging  effects  (in  the  form  of  line  aver- 
aging rather  than  volume  averaging)  have  been  investi- 
gated for  sonic  anemometers  (Silverman  1968,  Kaimal 
et  al.  1968,  Horst  1973a),  and  for  platinum  wire  ther- 
mometers (Kaimal  1968,  Wyngaard  1971).  They  have  also 
been  investigated  for  the  Lyman-Alpha  Hygrometer  by 
Andreas  (1981)  in  the  form  of  volume  averaging  over  a 
right  circular  cylinder. 

For  our  measurement  system  the  response  function 
associated  with  volume  averaging  (assuming  a  spherical 
volume)  is  given  by  the  following  equation,  which  is  de- 
rived in  appendix  C. 


10        a,  T*(k1R,x2,x3) 


[B3] 


where  R  is  the  radius  of  the  sphere  over  which  we  are 
averaging  and  is  the  wave  number  in  the  direction  of 
the  mean  horizontal  wind  and  by  Taylor's  hypothesis, 
k,  =  27rn/u,  where  u  is  the  mean  horizontal  wind  speed. 
Here  T»  (kjR.X2.X3)  *s  defined  as: 


T»(k1R,x2,x3)  = 


4x1 


J  lh  sin20  sin(k1Rx3cos^) 


ykjR  (l  +  x2Tj2  sin0)  dd 


[B4] 


where  JD  is  a  Bessel  function.  Equations  [B3]  and  [B4] 
were  integrated  numerically  (see  appendix  C  and  discus- 
sion below),  and  the  results  are  given  in  figure  B4,  which 
is  a  plot  of  Tv  (kjR).  This  figure  and  equation  [Bl]  show 
that  spatial  averaging  effects  are  important  only  when 
Tv  (k,R)  departs  significantly  from  a  value  of  unity, 
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Figure  B4.— Spectral  transfer  function,  Tv  (I^R),  associated  with 
spherical  volume  averaging  of  the  03,  N02,  and  NOx  concentra- 
tion fluctuation  measurements.  k,  =  2  71-n/u. 


which  occurs  whenever  k,  R  >  1  or,  equivalently, 
whenever  the  physical  dimensions  of  a  turbulent  at- 
mospheric eddy  are  comparable  to  or  smaller  than  the 
radius  of  the  averaging  volume.  Because  Tv  is,  in 
general,  a  function  of  wind  speed,  the  importance  of  the 
volume  averaging  effects  can  only  be  assessed  by  evalu- 
ating the  flux  loss,  equation  [Bl],  as  a  function  of  wind 
speed  with  and  without  Tv(n)  as  a  component  of  H  (n). 
This  will  be  discussed  below.  It  should  also  be  noted  that 
the  spherical  volume  over  which  we  are  averaging 
becomes  distorted  by  the  wind  speed;  however,  Taylor's 
hypothesis  allows  for  this  and  thus  k,  becomes  a  func- 
tion of  the  mean  wind  speed. 

Finally,  an  estimate  of  the  radius,  R,  is  needed.  A 
typical  value  of  R  can  be  estimated  from  measurements 
of  the  intake  flow  rate,  Q,  and  the  sampling  frequency, 
ns,  as  follows:  R  =  [3Q/(47rns)]1'3.  This  basically  assumes 
that  for  any  given  period  of  time,  At,  the  volume  QAt  is 
divided  into  nsAt  spherical  subvolumes  for  each  of 
which  we  make  a  concentration  measurement.  Given 
typical  flow  rates  of  about  150-170  1/min  and  a  sampling 
rate  of  10-20  Hz,  4  cm  is  a  reasonable  estimate  for  the 
radius  R.  This  value  will  be  used  in  the  next  section  for 
calculating  flux  loss  estimates. 


Gill  UVW  Anemometers 

Like  the  chemical  sensors,  the  3-axis  Gill  propeller 
anemometer  is  also  a  first  order  system;  hence,  the 
response  function  associated  with  this  instrument  is  also 
described  by  equation  [B2].  However,  the  response  time, 
t,  associated  with  the  Gill  propeller  is  not  a  constant  as 
it  is  with  the  chemical  sensors;  rather,  it  is  dependent 
largely  upon  the  horizontal  wind  speed,  u,  and  is  given 
by  t  =  D/u.  The  response  length,  D,  has  been  determined 
from  wind  tunnel  experiments  to  be  a  function  of  the 
attack  angle  of  the  wind.  For  a  propeller  measuring  the 
horizontal  wind  speed,  however,  this  effect  is  small,  and 
D  can  be  considered  constant  and  equal  to  1  m.  For  the 
vertical  propeller,  on  the  other  hand,  D  is  not  necessari- 
ly constant  but  will  increase  from  about  1  m  to  2  m  as 
the  flow  deviates  from  axial  to  about  80°  off  axis  and 
will  continue  to  increase  to  infinity  at  90°  off  axis.  There- 
fore, in  general  the  response  time  of  the  Gill  propellers 
varies  constantly  with  the  flow  velocity.  For  the  purposes 
of  this  work,  the  response  time  of  the  horizontal  pro- 
pellers is  assumed  to  be  1  m/u  and  that  of  the  vertical 
propeller  2.5  m/u.  Therefore,  the  response  function  for 
Gill  propellers,  which  uses  these  latter  two  response 
times,  is  an  approximation  to  the  true  response  function. 
For  more  details  on  the  subject  of  response  characteris- 
tics of  these  sensors,  see  Hicks  (1972)  and  Busch  et  al. 
(1980). 


Sensor  Mismatching 

In  addition  to  the  errors  in  flux  measurements  caused 
by  slow  response  times  of  individual  sensors,  mismatch- 
ing of  sensors  used  for  flux  measurements  can  introduce 
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further  errors.  Mismatching  errors  are  produced  when- 
ever sensors  have  different  response  characteristics. 
When  this  occurs,  there  is  a  relative  phase  shift  in  the 
time  series  gathered  by  the  different  sensors.  If  the 
relative  phase  shift  is  small,  then  the  loss  of  correlation 
(or  flux)  between  the  two  time  series  is  also  small. 
However,  as  the  relative  phase  shift  increases,  the  two 
time  series  can  become  increasingly  uncorrected  and, 
hence,  the  flux  measurements  may  be  seriously  under- 
estimated. In  general,  all  sensors  are  somewhat  mis- 
matched unless  they  happen  to  have  exactly  the  same 
response  characteristics  and,  therefore,  exactly  the  same 
response  functions.  However,  mismatching  is  most 
serious  when  flux  measurements  are  made  with  two  first 
order  sensors  that  have  very  different  response  times  or 
when  a  slow  response  sensor  is  used  in  conjunction  with 
a  fast  response  sensor.  In  our  case  this  corresponds  to 
(1)  the  chemical  flux  measurements  made  using  a  Gill 
propeller  during  low  or  high  wind  speeds  or  at  all  wind 
speeds  using  a  sonic  anemometer  and  (2)  the  heat  flux 
measurements  made  with  the  relatively  slow  responding 
vertical  Gill  propeller  and  the  much  faster  responding 
temperature  sensor. 

The  transfer  function  associated  with  mismatching  is 
given  simply  as  Tm(n)  =  cosf^-^J,  where  <p1  and  <p2  are 
the  phase  shifts  introduced  into  the  measurements  by 
sensor  number  1  and  number  2,  respectively.  For  first 
order  systems,  the  phase  shift  <p  =  tan-1  (27rn7),  from 
which  it  follows  that  cos<,c>  =  1/(1  +  (27rnT)2)1'2  and  sin^  = 
(27rnr)/(l  +  {2irnr)2)m.  Therefore,  using  these  latter  two 
relations  and  the  definition  of  Tm(n)  it  is  straightforward 
to  show  that: 


Tm(n)  = 


l  +  (2irn)2TjT2 


[(l  +  (27rnr1)2Kl  +  (27rnT2)2)]1 


[B5] 


Again,  this  transfer  function  is  appropriate  to  first  order 
systems  and  is  most  significant  whenever  r,  is  very  dif- 
ferent from  t2. 

Using  the  slower  responding  chemical  sensors  and  the 
faster  responding  vertical  sonic  anemometer  can  also  in- 
troduce mismatching  errors  in  the  03,  N02,  and  NOx 
flux  measurements.  For  this  case,  we  assume  that  [B5] 
is  valid  but  that  the  response  time  for  the  sonic  (e.g.,  tJ 
is  zero,  so  the  formulation  of  Tm(n)  simply  reduces  to 
1/(1  +  (27rnr)2)1'2,  where  r  is  the  response  time  of  the 
chemical  sensors.  This  approximation  will  produce  a 
slight  overestimation  of  the  flux  loss  because  the  sonic 
anemometer  does,  in  fact,  produce  a  small  phase  shift 
in  the  spectra  at  higher  wave  numbers. 


Sonic  Anemometers 


Kaimal  et  al.  (1968)  and  is  given  in  the  following 
equation: 


2r(4/3) 


r(i/2)  r(s/6)  (kjj2 


r 

J  o 


sin2  (k^x^) 


(x3/2); 


3/8 


(1+x2)  7/3       (l  +  x2)4/3 


dx3  [B6] 


where  lw  is  the  vertical  sonic  path  (=  25  cm  for  our 
sonic),  T  is  the  Gamma  function,  and  \  =  27rn/u  as  was 
discussed  above.  The  expression  given  in  [B6]  has  been 
adapted  from  Kaimal  et  al.  (1968)  to  permit  ease  in 
computation. 

The  response  function  associated  with  line  averaging 
of  the  horizontal  velocity  is  dependent  upon  the  instru- 
ment geometry  and  direction  of  the  horizontal  wind 
speed.  For  our  sonic,  which  uses  an  orthogonal  geometry 
in  the  horizontal  plane,  the  following  response  function 
will  be  used  to  describe  the  effects  of  line  averaging  of 
the  horizontal  velocity: 


Tl»  = 


sin2  (k1lu/2) 
(k1lu/2)2 


[B7] 


where  lu  is  the  horizontal  sonic  path  ( =  l^,  for  our  sonic) 
Figure  B5  shows  T\w  and  Tju  as  functions  of  kxl. 
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Figure  B5.— Spectral  transfer  functions,  T1  (k^)  and  T-,  (k^), 
associated  with  line  averaging  of  the  vertical  and  horizontal  wind 
velocity  fluctuation  measurements  by  the  sonic  anemometer. 

These  two  response  functions  assume  that  line  averag- 
ing takes  place  over  a  single  acoustic  path;  however,  our 
sonic  (and  many  others  as  well)  actually  employs  two 
parallel  acoustic  paths.  Horst  (1973b)  found  that  the  ef- 
fect of  path  separation  was  not  significant  in  the  vertical 
direction,  and  after  a  similar  investigation,  we  found  that 
for  our  sonic  the  path  separation  effects  were  negligible 
in  the  horizontal  as  well. 


Unlike  the  chemical  sensors  or  the  Gill  propellers,  the 
sonic  anemometer  is  not  limited  by  a  slow  response  time. 
Rather,  its  limitations  are  associated  with  line  averag- 
ing effects  (which  is  the  one-dimensional  equivalent  to 
the  volume  averaging  effects  of  the  chemical  sensors 
discussed  earlier).  The  response  function  associated  with 
line  averaging  of  the  vertical  velocity  is  taken  from 


Temperature  Sensor 

For  the  platinum  wire  resistance  thermometer  that  we 
are  using,  response  functions  for  both  line  averaging 
effects  and  the  response  of  a  first  order  system  will  be 
included. 


25 


The  response  time  of  this  thermometer  is  given  by 
Kaimal  (1968)  as  follows:  r  =  0.039/(1  +  V  3u),  where  u  is 
measured  in  meters  per  second.  Again  we  are  assuming 
that  this  sensor  is  a  first  order  sensor  and  we  use  equa- 
tion [B2]  for  describing  the  dynamic  frequency  response. 
The  response  function  appropriate  to  line  averaging  (of 
a  scalar  quantity  like  temperature  as  opposed  to  a  vector 
quantity  like  the  vertical  velocity)  was  first  developed  by 
Gurvich  (1962)  and  later  extended  by  Silverman  (1968). 
For  ease  of  computation,  their  results  have  been  adapted 
to  give  the  following  form: 


2r(4/3) 


r(i/2)  r(5/6)  (k,!)5 


sin2  (kjlx3/2) 


(x3/2): 


(l+x2)  4/3 


dx3  [B8] 


where  1  is  the  length  of  the  platinum  wire  (=  18  cm  in 
our  case).  Figure  B6  is  a  graph  of  T,  as  a  function  of  kd. 


Sensor  Separation 

For  a  variety  of  practical  reasons,  sensors  will  always 
be  separated  from  one  another.  Yet  this  separation  will 
also  cause  an  error  in  any  flux  measurement  for  which 
two  separate  quantities  are  being  correlated.  We  follow 
Moore's  (1986)  suggestion  that  the  following  response 
function  be  used  for  describing  both  lateral  and  longi- 
tudinal separation  effects: 
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Figure  B7.— Spectral  transfer  function,  Ts  (k.,s),  associated  with  the 
spatial  separation  of  any  two  eddy  correlation  instruments. 

Data  Logging  System 

Incorporated  into  our  data  logging  system  is  a  filter 
that  also  must  be  taken  into  account.  This  filter  is  a  high 
pass  recursive  filter  that  is  applied  to  all  raw  digitized 
sensor  data  to  remove  any  linear  or  low  frequency  long- 
term  trends  from  the  data.  These  long-term  trends  can 
adversely  affect  the  flux  measurements  (Kaimal  1975). 
The  transfer  function  associated  with  this  filter  is  ap- 
proximated by  Moore  (1986)  as  follows: 


trN  = 


(27rnTR)2a 
a  +  (2xnrR)2 


[B10] 


where  a  =  e-1'(T'ns)  and  r'  =  200  s  (for  our  filter),  and  tr 
=  a/[(l-a)nj. 


Ts(n)  = 


r(5/6) 


M5'6  K5/6(k,s) 


[B9] 


where  s  is  the  separation  distance  and  K5/6  is  the 
modified  Bessel  function  of  the  second  kind.  Figure  B7 
shows  a  plot  of  Ts  as  a  function  of  kjS.  As  Moore  (1986) 
cautions,  however,  [B9]  is  most  applicable  when  the 
separation  distances  are  not  too  great,  and  it  may  also 
lead  to  a  slight  overestimate  of  the  flux  loss. 
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Figure  86.— Spectral  transfer  function,  T1(k1l),  associated  with  line 
averaging  of  the  temperature  fluctuation  measurements  by  the 
platinum  wire  thermometer. 


Aliasing 

Aliasing  results  from  sampling  at  uniform  intervals  any 
continuous  signal  that  contains  energy  at  frequencies 
above  the  Nyquist  frequency.  The  Nyquist  frequency, 
which  is  defined  as  one-half  the  sampling  frequency,  will 
henceforth  be  denoted  by  n0.  As  a  result  of  aliasing, 
spectral  energy  at  frequencies  above  the  Nyquist  frequen- 
cy can  be  folded  back  into  the  spectra  at  frequencies 
below  the  Nyquist  frequency  at  harmonics  of  the  sam- 
pling frequency.  For  example,  the  effects  of  aliasing  upon 
measurements  of  Cwq(n)  can  be  written  as  follows 


(Kaimal  et  al 
H 


1968): 


(n)Cwq(n)  +  E 


JHwq 


0ns 


-n)cwq(jns- 


C~  in 

wq 


(n)  =  +  Hwq(jns  +  n)Cwq(jns  +  n)]  0  <  p  <  ^ 


[Bll] 


n>nr 


where  Cjjnjn)  is  the  measured  cospectra  of  w  and  q. 
Here  the  effects  of  the  sensor  response  functions  and  the 
filters  associated  with  the  data  logging  system  are  in- 
cluded as  components  of  Hwq(n).  The  summation  term 
of  [Bll]  is  the  aliased  portion  of  the  measured  cospec- 
tra. Although  not  specifically  included  in  [Bll],  any  spec- 
tra or  cospectra  can  also  be  contaminated  by  aliased  high 


26 


frequency  noise  associated  with  the  electronics,  nearby 
power  lines,  or  any  other  outside  power  source.  For 
sensors  with  short  response  times  (such  as  sonic  ane- 
mometers), antialiasing  filters  should  be  used  to  mini- 
mize the  aliasing  of  any  naturally  occurring  frequencies 
and  to  eliminate  any  high  frequency  contamination  of 
the  analog  signal  from  outside  sources.  These  antialias- 
ing filters  are  usually  designed  as  Butterworth  or  other 
similar  filters.  For  sensors  with  relatively  long  response 
times  (such  as  our  vertical  Gill  propellers  and  chemical 
sensors),  it  is  unnecessary  to  use  an  antialiasing  filter 
because  they  will  naturally  accomplish  these  goals. 

To  include  the  effects  of  aliasing  of  the  naturally  oc- 
curring frequencies  upon  the  flux  loss  computation,  it 
is  sufficient  for  our  purposes  to  confine  the  conse- 
quences of  aliasing  to  only  the  (n-ns)  term  of  [Bll]. 
Therefore,  we  can  simply  substitute  the  following  expres- 
sion for  the  integral  given  in  the  numerator  of  equation 
[Bl]: 

/o°Hwq(n)Cwq(n)dn  +  /^Hwqfrs-n^wqfrs-11^11 

However,  the  above  expression  can  be  further  simplified 
(by  a  change  of  variables)  to  the  following: 

/jHwq(n)Cwq(n)dn 

Therefore,  equation  [Bl]  can  now  be  written  in  the  fol- 
lowing form: 


AF 
F 


=  1- 


/^Hwq(n)Cwq(n)dn 


wq*. 


-CCwq(n)dn 


[B12] 


It  is  this  last  equation  that  was  used  to  estimate  the 
percentage  flux  loss  from  measurements  taken  with  our 
system.  The  sampling  frequency,  ns,  of  our  system  is 
about  14  Hz;  this  value  will  be  used  for  all  computations 
presented  herein. 

To  include  the  influence  of  the  antialiasing  filter  is 
quite  simple  with  [B12]  because  in  essence  these  filters 
are  designed  specifically  to  confine  the  effects  of  alias- 
ing to  only  the  (n-ns)  term  of  [Bll].  Therefore,  all  that 
needs  to  be  done  is  to  include  as  part  of  Hwq(n)  the  spec- 
tral response  function  of  a  particular  antialiasing  filter. 
For  example,  a  2-pole  Butterworth  filter  has  a  response 
function  given  as: 


TB(n)  =  l/(  !  + 


[B13] 


Other  types  of  antialiasing  filters  are  also  easily  included 
in  the  computations  for  the  flux  loss.  However,  before 
proceeding  to  the  flux  loss  computation,  we  need  to 
define  the  atmospheric  cospectra,  Cwq(n),  that  will  be 
used  in  [B12]. 

Atmospheric  Cospectra 

For  the  purposes  of  this  work,  the  cospectra  of  momen- 
tum and  sensible  heat  are  the  same  for  stable  atmos- 
pheric conditions  and  are  given  by: 


Cuw(n)  =  CwT(n)  = 


1  +  2.34(z/u)21n2 


[B14] 


where  z  is  the  height  of  the  sensors  array  above  the  zero 
plane  displacement  of  the  surface  (z  =  6  m  for  our  site). 
This  expression  is  taken  from  Kaimal  et  al.  (1972),  and 
although  the  normalization  is  not  necessarily  standard, 
for  our  purposes  nothing  more  is  needed.  For  unstable 
atmospheric  conditions,  the  cospectra  of  momentum  and 
sensible  heat  are  different  and  given  as  follows: 


Cuw(n) 


20.78 


= 


CwT(n)=  , 


[l  +  31(z/u)n]1575 

12.66 
[l  +  9.6(z/u)n]24 

12.92 
[l  +  26.7(z/u)n]1375 
4.378 
[l  +  3.8(z/u)n]24 


(z/u)n  <  0.24 


(z/u)n  >  0.24 


(z/u)n  <  0.54 


(z/u)n  >  0.54 


[B15] 


[B16] 


These  latter  two  equations  were  taken  from  Moore 
(1986). 

For  correcting  the  dry  deposition  fluxes  of  03,  N02, 
and  NOx  ,  we  will  assume  that  as  atmospheric  scalars, 
their  cospectra  are  the  same  as  that  of  temperature. 
There  is  some  experimental  evidence  for  the  similarity 
between  the  spectra  and  cospectra  of  temperature  and 
the  spectra  and  cospectra  of  other  scalars  such  as 
humidity  and  C02  (Redford  et  al.  1980,  Smith  and 
Anderson  1984,  Ohtaki  1985,  Anderson  et  al.  1986). 
However,  no  direct  measurements  of  the  spectra  or 
cospectra  of  O,,  NO, 


or  NOx  are  available. 


Atmospheric  Spectra 

With  two  simple  modifications  to  [B12],  the  variance 
losses  due  to  sensor  limitations  can  be  estimated  in  a 
manner  completely  analogous  to  that  used  for  estimating 
the  flux  losses.  First  Hwq(n)  is  replaced  by  the  ap- 
propriate transfer  function  for  a  given  instrument  (i.e., 
Hw  or  Hc  or  HT  replaces  Hwq)  and,  second,  Cwq(n)  is 
replaced  by  the  spectra  appropriate  to  the  quantity  be- 
ing measured. 

For  stable  atmospheric  conditions,  the  following 
variance  spectra,  taken  from  Moore  (1986),  will  be  used: 


SqqN  = 


l  +  aQ[(z/u)n]5'3 


[B17] 


where  q  refers  to  vertical  velocity  or  temperature  or  any 
trace  gas  (03,  NOz,  or  NOx).  Again,  the  spectra  asso- 
ciated with  the  trace  gases  is  assumed  to  be  the  same 
as  that  associated  with  temperature.  For  q  =  vertical 
velocity,  aq  =  3.124[0.838  +  1.172(z/L)]-5'3,  and  for 
q  =  temperature,  aq  =  3.124(0.0961  +  0.644(z/L)3'5]-5'3.  Here 
z/L  is  the  measure  of  atmospheric  stability  and  is  defined 
as  -kgzw'0'/(u;|T)  where  L  refers  to  the  Monin-Obukhov 
length,  k  is  the  von  Karman  constant  (  =  0.41),  g  is  the 
acceleration  due  to  gravity  (  =  9.8  m/sec2),  w'6'  =  1.05 
w'T'  at  our  site  with  6  being  the  potential  temperature, 
and  u»  is  the  friction  velocity  (u,  =  (-u'w')1'2).  For  stable 
conditions  z/L  is  positive,  and  for  unstable  conditions 
z/L  is  negative.  Neutral  stability  refers  to  z/L  =  0. 
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For  an  unstable  or  neutral  atmosphere  (z/L  <  0),  the  ver- 
tical velocity  variance  spectra  is  taken  from  H6jstrup 
(1981)  and  the  temperature  variance  spectra  is  taken  from 
Kaimal  et  al.  (1972). 


,(n)  = 


STT(n)  = 


l  +  5.3[(z/u)n5/3] 

14.94 
[1  +  24(z/u)n]5/3 

6.827 
[1  +  12.5(z/u)n]5'3 


16(-z/L)2/3 
[l  +  17(z/u)n]5'3 

(z/u)n<0.15 
(z/u)n>0.15 


[B18] 


[B19] 


With  all  the  above  simplifying  assumptions  in  mind, 
we  applied  these  response  functions  and  atmospheric 
cospectra  to  our  eddy  correlation  system.  For  the  vertical 
Gill  propeller,  Hw(n)  =  T^(n)TR(n)  where  the  w  subscript 
and  superscript  refer  to  the  vertical  velocity  sensor  and 
Hw(n)  is  the  spectral  response  function  for  this  instru- 
ment. Likewise,  for  the  horizontal  propellers,  Hu(n)  = 
Tj(n)TR(n);  for  the  temperature  sensor  HT(n)  =  Tj(n) 
Ti(n)TB(n)TR(n),  and  finally  for  the  chemical  sensors 
HJn)  =  T§(n)TR(n)Tv(n).  For  the  vertical  sonic  Hw(n) 
=  Tiw(n)TB(n)TR(n)  and  for  the  horizontal  sonic  Hu(n) 
=  Tiu(n)TB(n)TR(n).  Therefore,  the  total  system  response 
functions  for  heat  and  momentum  and  chemical  fluxes 
are  given,  respectively,  as  follows: 

HwT(n)=  V  Hw(n)HT(n)  T™T  (n)  T£T(n) 

Huw(n)=  V  Hu(n)Hw(n)  T»w  (n)  T»w(n)       j  [B20] 

Hwc(n)  =  \f  Hw(n)Hc(n)  T™c  (n)  T™c(n) 

where  the  instrument  separation  and  mismatching  ef- 
fects are  denoted  by  Ts(n)  and  Tm(n),  respectively,  and 
the  superscripts  refer  to  that  particular  pair  of  instru- 
ments. It  is  important  to  note  here  that  [B20]  is  an  approx- 
imation. With  [B20]  we  have  "constructed"  a  response 
function.  However,  it  is  more  appropriate  to  derive  a 
response  function  for  correcting  the  flux  measurements 
similar  to  the  way  [B3],  for  example,  was  derived  for  cor- 
recting variance  measurements;  but,  this  more  funda- 
mental approach,  except  for  the  preliminary  efforts  of 
Kristensen  and  Fitzjarrald  (1984),  remains  largely 
undeveloped  both  theoretically  and  observationally. 


Computational  Methods 

Several  different  numerical  methods  were  employed 
to  evaluate  equation  [B12].  The  denominator  of  this  ex- 
pression was  evaluated  analytically  for  the  cospectra 
given  by  [B14],  [B15],  and  [B16].  The  numerator  was 
evaluated  numerically  using  the  Romberg  method  of  inte- 
gration (Dahlquist  and  Bjorck  1974).  The  different  trans- 
fer functions  used  for  Hwq(n)  in  [B20]  were  represented 
either  directly  by  a  specific  analytical  expression  (e.g., 
[B2],  [B5],  [BIO],  and  [B13])  or  indirectly  by  a  numerical 
table  (e.g.,  [B3],  [B6],  [B7],  [B8],  and  [B9]).  The  numerical 


tables  were  provided  by  a  simple  evaluation  of  a  func- 
tion (e.g.,  [B7]  and  [B9]),  or  by  the  numerical  evaluation 
of  an  integral  (e.g.,  [B3],  [B6],  and  [B8]).  As  with  [B20], 
these  latter  three  integrals  were  also  evaluated  with  the 
Romberg  method  of  integration.  Finally,  the  numerical 
evaluation  of  [B20]  sometimes  required  interpolating  bet- 
ween points  of  the  numerical  tables  representing  [B3] 
and  [B6]  through  [B8].  This  was  done  using  a  quartic 
polynomial  interpolator. 


APPENDIX  C. 
Transfer  Function  for  Volume  Averaging 

By  using  an  intake  tube  to  draw  atmospheric  samples 
into  a  reaction  chamber,  our  chemical  sensors  average 
the  turbulent  fluctuations  of  03,  N02,  and  NOx  over  a 
spherical  volume.  Therefore,  the  spectral  response  of 
these  sensors  must  degrade  at  higher  wave  numbers. 
This  appendix  provides  a  derivation  of  the  response 
function  associated  with  the  volume  averaging  effects 
of  the  intake  tube  (Tv(n)  in  equation  [B3]),  which  can  be 
used  if  necessary  to  correct  the  flux  measurements  for 
these  effects.  We  begin  by  considering  only  spatial 
averaging  effects;  therefore,  we  will  ignore  other  limita- 
tions of  the  sensors  (such  as  a  relatively  long  response 
time),  which  also  influence  the  spectral  response  of  these 
instruments. 

From  Andreas  (1981)  it  is  straightforward  to  show  that 
the  following  relation  holds  between  the_  measured 
constituent  (03,_N02,  or  NOx)  spectra,  <£m(k),  and  the 
true  spectra,  </>(k): 


</>m(k)  = 


(4ttR3)2 
3 


vol 


ik.s. 


'lift 


-ik.s 


ds' 


m  [ci] 


where  k  is  the  three-dimensional  wave  number  and  R 
is  the  radius  of  the  volume  over  which  we  are  averag- 
ing. In  arriving  at  [Cl]  Andreas'  (1981)  results  for  a  right 
circular  cylinder  have  been  modified  to  an  expression 
appropriate  to  a  sphere.  We  have  chosen  the  same  right- 
handed  coordinate  system  as  he  did  (x  =  direction  of  the 
mean  horizontal  wind;  y  =  direction  in  the  horizontal 
plane  perpendicular  to  the  x  axis;  z  =  vertical  direction 
=  alignment  direction  of  the  intake  tube).  However,  as 
in  his  case,  the  choice  of  this  particular  coordinate 
system  is  immaterial. 
From  [Cl]  it  is  clear  the  integral,  I  =  Je  lksds£  must 

be  evaluated.  In  cylindrical  coordinates  tne  integral  is 
given  by  the  following  expression: 


1=  [27rfRf+,!Jr  ^),!,'2eik3Zeir{kicose  +  k2sin0d0rdrdz  [C2] 
Jo  Jo  J -(R-z2)1'2 

where  klf  k2,  and  k3  are  the  three  components  of  k  along 
the  x,  y,  and  z  axes.  The  integral  over  6  can  be  done  ex- 
plicitly (Andreas  1981)  as  can  the  integral  over  z  (Silver- 
man 1968,  Wyngaard  1971),  which  yields: 


:  =  2tt 

J  n 


sin((R2-r2)1/2  k3) 
V2 


J0(r(k2  +  k2nrdr  [C3] 
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where  J0  is  a  Bessel  function. 
Then  changing  variables  from  r  to  -  =  sin  6  yields: 


1=  2»R»  J 


(k3R)  o 


sin20  sin(k3  Rcos0)  Jo(Rsin0(k2  +  k22yl2)dd 


Finally,  using  [Cl]  and  [C4],  we  have 
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T'  (k,R)  = 


4(k3R)2 


V/2 

sin20  sin(k3Rcos0) 
°    J0  (Rsin0(k2  +  k22)m)  dd 


[C4] 


[C5] 


where  T'  (k,R)  is  defined  as  <£m(k)/<£(k). 

However,  only  the  one-dimensional  spectra  is  observ- 
able and  the  true  and  measured  one-dimensional  spec- 
tra (denoted  by  ^(kj  and  ^rm(k1),  respectively)  are 
related  to  </>(k)  as  follows: 

*(kj  =  ft  4>(k)dk2dk3  [C6] 

 CO 

*miK)=  FS  </>m(k)dk2  dk3=  Ff  T'(k,R)0(k)dk2dk3  [C7] 


The  transfer  function  we  seek,  Tv(n),  is  defined  as 
*m(k1)/1Jr(k1).  With  this  in  mind,  the  integrals  for  Tv(n) 


(viz.,  C5,  C6,  and  C7)  can  be  computed  once  cj>(k)  is 
specified.  Assuming  local  isotropic  turbulence,  the  at- 
mospheric variance  spectra  for  a  scalar  quantity  is  given 
simply  as  0(k)  =  Ak  n'3  where  A  is  a  constant  and  k  is 
the  magnitude  of  k.  Therefore,  Tv(n)  can  be  defined  as 
follows: 


Tv(n)  = 


F[  k-»'3T'  (k,R)dk2  dk3 

 CO 

ft  k-11/3dk2dk3 


[C8] 


The  denominator  can  be  evaluated  analytically  to  give 
6tt/5  k-5/3.  With  this  last  result  [C8]  and  [C5]  are 
equivalent  to  equations  [B3]  and  [B4]  of  appendix  B  and 
all  that  remains  is  to  note  that  x2  =  k2/k1  and  x3  =  k3/kj 
and  to  change  variables  of  integration  appropriately. 

The  integrals  given  in  [B3]  and  [B4]  were  evaluated 
numerically  using  the  Romberg  method  of  integration 
(Dahlquist  and  Bjorck  1974).  The  Bessel  function  J0was 
also  evaluated  numerically  using  rational  functions  and 
rational  polynomial  approximations  of  J0  (Press  et  al. 
1986).  Finally,  because  evaluating  Tv(n)  required 
evaluating  a  triple  integral,  some  care  in  design  and 
testing  of  the  computer  code  was  necessary  in  order  to 
assure  that  figure  B4  is  a  good  representation  of  Tv(n). 


APPENDIX  D.  Table  of  Symbols 


Symbol  Definition  Unit 

c  pollutant  concentration  kg/m3  or  ppb/v 

Cwq(n)  atmospheric  cospectrum  of  w  and  q  as  a  function  of  n 

D  Gill  propeller  anemometer  response  length  m 

F  corrected  pollutant  flux  kg/m2/s  or  (ppb/v)(m/s) 

Fm  measured  pollutant  flux  kg/m2/s  or  (ppb/v)(m/s) 

g  acceleration  due  to  gravity  m/s2 

Hwq(n)  total  cospectral  response  function  for  vertical  wind  w  and  quantity  q  as  a 
function  of  frequency 

k  von  Karman  constant 

kj  spectral  wave  number  in  ith  direction  cycles/m 

L  Monin-Obukov  length  m 

n  frequency  s'1 

q  specific  density  of  an  arbitrary  measured  quantity 

Q  flow  rate  liters/min 

rs  surface  resistance  s/m 

rt  total  transfer  resistance  s/m 

r..„  correlation  coefficient  between  w  and  c 
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R  radius  m 

Sqq(n)  variance  spectra  for  any  quantity  q  as  a  function  of  n 

t  time  s 

T  temperature  °C 


TB(n)  response  function  for  Butterworth  filter 

Td(n)  response  function  for  dynamic  frequency  response 

Tj(n)  response  function  for  line  averaging 

Tm(n)  response  function  for  mismatched  sensors 

TR(n)  response  function  for  high-pass  recursive  filter 

Ts(n)  response  function  for  sensor  separation 

Tv(n)  response  function  for  volume  averaging 


u  horizontal  wind  speed  m/s 

u*  friction  velocity  m/s 

vd  deposition  velocity  cm/s 

w  vertical  wind  speed  m/s 

x;  kj/kj,  i  =  2,3 

z  vertical  coordinate 

z/L  stability  parameter:  measurement  height  divided  by  Monin-Obukov  length 

z0  roughness  length  cm 

Oj  temperature  variance  (°K)2 

ct£  variance  of  c  (ppb/v)2 

(j^  vertical  velocity  variance  (m/s)2 

0  potential  temperature  °K 

pc  density  of  pollutant  c  g/m3 

t  sensor  response  time  s 

(p  phase  angle  due  to  sensor  mismatch 

\pc  universal  stability  influence  function 
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ical  theory;  and  (c)  micrometeorological  site  characteristics. 

Keywords:  Eddy  correlation,  dry  deposition,  Pawnee  Grasslands,  03, 
N02,  NOx,  micrometeorology 
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U.S.  Department  of  Agriculture 
Forest  Service 

Rocky  Mountain  Forest  and 
Range  Experiment  Station 


The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 


RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 


Albuquerque,  New  Mexico 

Flagstaff,  Arizona 

Fort  Collins,  Colorado* 

Laramie,  Wyoming 

Lincoln,  Nebraska 

Rapid  City,  South  Dakota 

Tempe,  Arizona 


'Station  Headquarters:  240  W.  Prospect  St.,  Fort  Collins,  CO  80526 


